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Abstract: We study the evolution of an energetic jet which radiates gluons while propa- 
gating through a dense QCD medium modeled as a random distribution of color sources. 
Motivated by the heavy ion experimental program at the LHC, we focus on the medium- 
induced radiation of (relatively) soft gluons, which are abundantly emitted at large angles 
and thus can transport a small fraction of the jet energy far away from the jet axis. We 
perform a complete calculation of the medium-induced gluon branching in the regime 
where the gluons that take part in the branching undergo multiple soft scattering with the 
medium. We extend the BDMPSZ theory of radiative energy loss by including the trans- 
verse momentum dependence in the kernel that describes the branching and by analyzing 
the correlations between the two offspring gluons. We demonstrate that these gluons lose 
color coherence with respect to each other over a time scale that is comparable to the 
duration of the branching. It follows that interference effects between successive emissions 
are suppressed, a necessary ingredient for a description of multiple emission of soft gluons 
by a probabilistic, branching process. 
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1 Introduction 

The LHC heavy ion program has produced a wealth of remarkable results [1-4] that mo- 
tivate new theoretical efforts in the study of jet propagation in matter. For one thing, 
these results confirm those obtained at RHIC [5-7], indicating that matter formed in the 
collisions strongly suppresses the yield of high-p^ hadrons as compared to the yield that 
one deduces from proton-proton collisions after proper scaling by the number of binary 
collisions in the nucleus-nucleus collision. This suppression of high-pr hadrons, referred to 
as "jet-quenching", is usually attributed to the energy loss of the leading partons caused 
by the radiation of soft gluons induced by their collisions with the matter constituents. 
The theory of radiative parton energy loss has been developed in late 90 's. It is commonly 
referred to as BDMPSZ theory, from the names of the original authors [8-12]. Further 
developments are presented in [13-20]. In this theory, the energy loss is characterized by 
a single parameter, a transport coefficient called q (the 'jet quenching parameter'), which 
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measures how much transverse momentum AA;^ a given parton acquires through multiple 
scattering as it travels through the medium over a distance A/: A/c^ = q/S.1. 

However, the LHC experiments provide much more detailed information about what is 
going on as the jet propagates through the medium, beyond the mere evidence for energy 
loss. There is in particular clear evidence that the jet shape is affected with a large number 
of soft particles being emitted at larger angles, i.e., outside the jet cone, and carrying a 
small (some ten percent or so) of the total energy [2~4]. While these soft particles do not 
contribute much to the energy loss, they carry important information, as we shall see, on 
the basic microscopic mechanisms at work. In order to fully exploit this new available 
information, a more complete description of exclusive jet observables (such as jet-shapes, 
particle correlations, etc) is called for. This paper presents the first step towards such a 
more complete theory of jet propagation in matter, which is valid in a specific regime that 
we shall specify shortly. 

In order to put our work into perspective, we need first to recall a few basic features of 
the BDMPSZ theory [8~12]. When propagating through a medium in which it undergoes 
multiple scattering, a high energy parton can radiate a gluon over a typical time scale 
("formation time") given by 

-1 1-1 

where cj and k±_ are respectively the energy of the radiated gluon and its transverse mo- 
mentum (with respect to the parent gluon). The quantity 1/tj, which, as suggested by 
Eq. (1.1) may be read as a non relativistic energy, with a; playing the role of a mass (the 
relevance of this analogy will be clarified in the main text) is essentially the amount of 
energy that is required to put the gluon on-shell. In the medium, such an energy is ob- 
tained from multiple scattering with the plasma constituents, each collision providing the 
colliding hard partons some transverse momentum. As mentioned earlier, the rate at which 
transverse momentum is accumulated by a parton along its trajectory is given by that is 
A/c^ = q/S.t. One sees therefore that emission of a gluon of a given energy u) can take place 
if the transverse momentum acquired during matches 2uj/t^, that is if l/r^ ~ [qrjuj). 
This provides a self-consistency condition that determines, as a function of w, the time 
scale for in- medium splitting, that we shall denote by r^^,: 




~ \Hr- (1-2) 



Related to r^^ , it is also convenient to define k^^ , the typical transverse momentum acquired 
during the time r^^.: fc^^ = qr^^^ or k^^(oj) ~ (2a;g)-^/^. Note that the time scale of the 
branching process is larger for harder gluons, because it requires more collisions to put a 
hard radiated gluon on-shell than a soft one {k^^ is as slowly growing function of w). Note 
also that the radiation is emitted with a characteristic angle 0^^ ~ ^hJ'^ ~ (^Z'^'^)^^^- Thus 
the softer the emission, the larger the emission angle, with all angles such that 6 > 9c, the 
minimal angle 6^, = {qL^) corresponding to r^^ = L, the length of the medium. 
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The BDMPSZ energy spectrum of the radiated gluons is of the form 



w— ~ W— = aW— , 1.3 

da; vr V V 

for uj < ujc (and more strongly suppressed for uj > uJc)- The frequency Dc is that for 
which T^^^{ujc) ~ L, that is, Uc = ^qL"^- The first factor in Eq. (1.3) is the standard 
bremsstrahlung spectrum for radiation by the parent gluon. The correction factor, which 
we can write as s/uiiJZ: = L/t^^^{ijj), increases as ui decreases, and the spectrum (1.3) 
should be cut-off at a minimal frequency o^bh' which multiple scattering ceases to be 
important and the radiation is produced by incoherent collisions (Bethe-Heitler spectrum) : 
iOg^ is the frequency for which the formation time is of the order of the mean free path 
i between successive collisions, that is ~ ^- Thus, L/t^^{uj) is the number of 

effective scattering centers, which is maximum for uj ~ uj^^ and is of order 1 in the vicinity 
of cjJc- The decrease of L/t^_^{ijj) as u increases may be understood as a consequence of the 
Landau-Pomeranchuk-Migdal (LPM) effect. 

This mechanism for gluon production, and the approximations used to calculate it, 
require the formation time to be much larger than the mean free path but smaller than 
the size L of the medium, that is oj^^ <^ u < ujc- For energies within this range, there is a 
large number of scattering centers, of order t^^(u})/£ ^> 1, which coherently contribute to 
the emission process. This can be rephrased in terms of the typical transverse momentum 
exchanged in one collision. Having in mind a picture of the medium where the typical 
collisions are induced by a screened one gluon exchange, we call that typical momentum 
niD (with niD the screening mass). Then, from the definition of q given earlier, rnjj = q£, 
and the condition oj^^ <^ u) translates into k'^^i^j) ^» mj^. 

The initial motivation for the BDMPSZ theory was to provide a framework for calcu- 
lating the energy loss. An estimate of this energy loss can be obtained by integrating the 
spectrum (1.3): 

AE = / da;a;— — ~ aujc ~ aqL^. (1-4) 

Ju^o da; 

This is dominated by the upper limit a; = ojc, the maximal energy that can be taken 
away by a single gluon. Such emissions of hard gluons are rather rare: they occur with a 
probability of order a. However, from the spectrum (1.3) one can also compute the average 
number of gluons emitted with energies larger than a given value oj: 

. , , . ,dN L , , 

AN(u}) = da;' — ~ a—-. 1.5 
da;' Tf{uj) 

As long as AA^(a;) < 1 (that is, as long as a; > uic), it may be identified to the probability 
to emit one gluon with energy a;' > a;. In such a case, the probability for multiple emissions 
is small. This is the case of the relatively hard emissions that dominate the energy loss. 
But for sufficiently soft gluon emissions, such that 

a; < a;s = a;c , (1.6) 
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one has AA^(a;) > 1 and then the multiple emissions are clearly important. Note that these 
multiple emissions are 'large angle' emissions, with 9 ^ a course this can 

only occur if the medium is large enough, since the condition co^^ <C ujs implies i <C aL. 

To summarize, there is a regime, characterized by gluon energies in the range uj^^ ^ 
oj <ti iOc, were medium induced radiation dominates, and where multiple emissions are 
important. This is the regime that we explore in this paper. Note that the conditions on u 
and L that characterize this regime are not very restrictive, and in fact they are expected 
to be well satisfied in the LHC experiments. For a rough orientation, taking q ~ IGeV^/fm 
and L ~ 6 fm, one finds Uc ~ 100 GeV, while the typical energy of soft particles is in the 
GeV range. 

When multiple emissions become important, i.e., when AA^(a;) > 1, the whole calcula- 
tionnal setting needs to be revised. In technical terms, when aL/r^^ ~ 1 , the perturbative 
expansion breaks down, and powers of aL/r^^ have to be resummed. (For instance, the 
longitudinal phase-space for two independent successive branchings goes like and the 
process is of order [aL / t^_^)'^ .) The treatment of multiple emission is a priori complicated 
by interferences between various high order processes. Various aspects of interference phe- 
nomena for medium-induced gluon radiation have been recently studied, but only for the 
case of a frozen configuration of the emitters (a pair of partons forming a 'colour antenna') 
[21-25]. However, in the soft regime of interest for us here, these interferences turn out 
to be negligible. Indeed, as we shall see, color correlations between the offspring gluons 
disappear on the same time scale r^^ at which the splitting occurs. As a consequence, the 
newly formed gluons propagate independently from each other except for a relatively small 
period r^^ ^ L. Subsequent emissions from these gluons can interfere with each other 
only if they occur during that period of order r^^ where color coherence is still present. 
This implies that the longitudinal phase space for interference effects is smaller, by a factor 
Ty^^/L <^ 1, than the corresponding phase space for independent emissions. Note that these 
features are specific to medium induced radiation^. By contrast, in vacuum, color is con- 
served along the parton shower, which implies coherence of successive parton branchings, 
that is responsible for angular ordering [26-29]. 

The argument of the previous paragraph suggests that multiple emissions could be 
treated as a probabilistic cascade of independent branchings, since this naturally resums 
the terms with maximal powers of [aL/r^^). There is however a limit to this argument, since 
successive emissions may overlap with each other when the energies of the produced gluons 
become too soft. Consider indeed the probability for emitting one gluon with frequency 
uj > during an interval At. In the regime where this probability is small, this may be 
estimated from the BDMPSZ spectrum, Eq. (1.3), and is (see Eq. (1.5)) 

P(At;L.o) ~ a^^. (1.7) 

^In fact, the same estimate for the color decoherence time, that is, a time scale of order t^^^{lo), would 
also follow from the previous analyses of a quark-antiquark color antenna in Refs. [21, 23] provided in 
these analyses one identifies the angular opening of the antenna with the formation angle 9^^^ ~ [q/ uf^^^'^ ■ 
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This becomes of order unity when At = rrad(wo) with 




a 



1 




(1.8) 



This quantity Ti.ad('Wo) rnay be understood as the typical time interval that a gluon with 
energy u; can survive without emitting any gluon with energy uj' within the interval a;o < 
w' < w. Note that it is the emission of the softest allowed gluons which controls the size 
of this interval. Because of the inverse power of a in Eq. (1.8), ri.ad('^o) is parametrically 
larger than the formation time of the radiated gluon. This property would hold all along the 
cascade if produced gluons are emitted with comparable frequencies: successive emissions 
remain well separated, the phase-space for overlapping emissions is small and successive 
emissions can be treated as independent. This is what is implicitly assumed in the previous 
paragraph. Under these conditions, the whole branching process reduces to a classical 
stochastic process^ obtained by iterating the elementary building block corresponding to 
one single splitting of a gluon into two gluons. This is the process that will be studied in 
detail in this paper. 

However, one cannot exclude a priori the situation where a very soft gluon is emitted, 
with Ti.ad('Wo) comparable to the formation time of the parent gluon (with energy uS). The 
condition Ti.ad(wo) > ''"br('^) holds if ^ o^uj. Thus, if cjq is parametrically smaller than 
uj (by two powers of o;), the interval Tj-ad (wo) between two successive emissions becomes 
of the order of the duration of the branching process that has produced the gluon uj in 
the previous step, and then the argument of independent emissions breaks down. This 
particular issue related to the role of very soft emissions in an actual cascade will be 
discussed in a forthcoming publication. 

The calculation to be presented in this paper is technically involved, but conceptually 
simple. We consider a high energy gluon propagating in a quark- gluon plasma. The 
medium is modeled by a random color field whose fluctuations account for collisions with 
plasma constituents and determine q. In Section 2, we present the general strategy of the 
calculation by discussing the well known phenomenon of momentum broadening. Then, in 
Section 3 we discuss the general structure of the gluon splitting amplitude and probability. 
The proof of factorization of the splitting probability requires a special study of a 4-point 
function, that we shall be able to explicitly compute only in the limit of a large number 
of colors A'c ^ 1. The corresponding analysis is presented in Section 4. However, the 
physical argument for the color decoherence between the offspring gluons is quite general 
and we expect it to remain valid for any value of ISi^- Section 5 is devoted to the calculation 
of the splitting kernel and the completion of the formula for the splitting probability. In 
the conclusion we summarize the results and their range of validity. Technical material is 
gathered in the Appendices. Appendix A reviews general features of gluon propagation 
in a specific background field. Appendix B is devoted to the calculation of specific path 

^This conclusion brings some support to previous phenomenological studies inspired by the BDMPSZ 
theory [30-32] and in particular the MonteCarlo codes developed in Refs. [33-35], which have already 
assumed the factorization of subsequent emissions. 
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integrals that enter the n-point functions that are studied in the text. Finally Appendix 
C provides details on the color algebra that is needed to calculate the 4-point function. 

2 Setting up the calculation with a simple example 

We shall set up formalism, and fix the notation, by reviewing the well known phenomenon 
of momentum broadening of a parton traveling through a quark-gluon plasma. This will 
also give us the opportunity to present our main results concerning the medium-induced 
gluon branching, anticipating on the detailed calculations to be done in the forthcoming 
sections. 

We consider an energetic parton propagating through a quark-gluon plasma. For 
simplicity, we restrict ourselves to the case where this energetic parton is a gluon. The 
extension to the case where it is a quark is straightforward. The gluon is produced inside 
the medium, via some hard scattering process, and then it propagates along a distance L 
until it escapes into the vacuum. During this propagation, the gluon interacts with the 
medium, and exchanges with it color and transverse momentum. 

We assume that the energetic parton propagates at nearly the speed of light along the 
axis and we use light-cone coordinates, e.g. = {x~^,x~,x), with'^ 



together with the light-cone gauge A'^ = 0. We describe the medium as a random color 
field with correlation function 



where the angular brackets denote the medium average. More general situations may 
be considered, e.g. we may allow 7 to depend on x~^, but we shall not consider such 
straightforward extensions here. 

In writing Eq. (2.2), we have made implicitly several simplifications. The coupling of 
the energetic gluon with the medium is described in the eikonal approximation, i.e., we 
assume that it couples only to the component A~ of the gauge field. Furthermore, this field 
is probed only at small values of x~ (closed to the trajectory of the gluon x~ ^0), and we 
can ignore its dependence on x~ . In other words, we assume that the + component of the 
gluon momentum is conserved during the propagation through the medium. Furthermore, 
because of Lorentz time dilation, the hard gluon has a very poor resolution in x'^, and the 
medium correlations, that have a finite extent in x^, appear to it as eff'ectively local : this 
is the origin of the 5-function 6{x~^ — y~^) in the correlator (2.2). By the same token, the 
field of the medium, which has a finite longitudinal extent L, appears frozen during the 
time when the hard gluon traverses it. That implies that the average over the medium will 
be done after squaring the amplitudes. We take this average over the fields to be Gaussian, 

^From now on, we renounce to the subscript _L on transverse components, that is, we write e.g. x± = x 
to alleviate notations. 
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assuming that the corrections to this approximation are of higher order in the coupHng 
strength. Finally, as obvious in Eq. (2.2) we assume homogeneity in the transverse plane, 
i.e., the correlation 7 is a function of x — y alone. In summary, the problem that we 
address is the propagation of an energetic gluon in a random background A~ field that 
is independent of x~ , with Gaussian correlations. Technical aspects of this problem are 
studied in detail in Appendix A. 

Although this will not enter explicitly in our derivations, it is perhaps useful to have in 
mind a specific model for the medium. We then briefly discuss the case where this medium 
is a weakly coupled quark-gluon plasma in thermal equilibrium at high enough temperature 
T. In this case, the medium constituents are quarks and gluons with energies and momenta 
p ^ T. Assuming that the charge carriers are correlated over distances determined by 
the screening length (inverse Debye mass) rn~^, one can estimate the correlator of field 
fluctuations: 

jix -y)=9'n[^^ e^'^<-y) (2.3) 



(27r)2 {q^+m 



D ' 



where re oc is the density of (point-like) color charges (weighted with appropriate color 
factors), and 1/(q^ + ™'^) is the screened Coulomb propagator. In fact, the quantity which 
will enter our analysis is not 7 itself, but the difference 7(0) — 7(r), 

7(0) - 7(r) =g'n / ~ 'T,2 • (2-4) 

J (27r)^ (q^ + rei^)^ 

The integral over q in Eq. (2.4) is dominated by g < 1/r, and in the relevant case where, 
typically, 1/r ^ m^, it is logarithmically sensitive to all the scales within the range 
rei,^ <C g ^ 1/r. To leading logarithmic accuracy, it can be evaluated by expanding the 
complex exponential to second order. One gets 

7 0) -7(r ~ -fir rer -r^^ 2x2^ ^9 nr In^-^. (2.5 

In the context of the present calculation the relevant values of r are typically fixed by 
transverse momentum broadening. In the regime dominated by multiple scattering, one 
often uses the harmonic approximation, where one ignores the weak dependence on r of 
the logarithm, and set 

g^N,[^{0)-j{r)]^^qr', (2.6) 
with the jet quenching parameter q (evaluated at the typical scale 1/f) given by 

q ^ ATralNcTiln^^ . (2.7) 

Parametrically, we have, ~ gT, i ~ l/g'^T, q ~ g^T^, and u)bh ~ T. 

We now return to our main discussion. We assume that the energetic gluon is created 
inside the medium at time x"*" = to by some local current J (to), and then propagates up 
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to time x"*" = tL, when it escapes the medium. Note that, as mentioned akeady in the 
introduction, we ignore in our calculation the 'vacuum radiation' that would be associated 
with the possible virtuality of the energetic gluon. (In what follows we shall often denoted 
the light-cone time simply as t, to alleviate the notation. Note that tL = V^L.) The 
probability amplitude for finding the (on-shell) gluon at t^, with momentum = {k~ = 
fc^/2fe+, k~^, fc), colour b, and polarization A, can be obtained from a standard calculation, 
as explained in Appendix A (see also Refs. [21, 23, 36] for previous applications of a similar 
formalism). One gets 



M\{k+,k) 



(2vr)2 



{k b\ g{tL, to) \po a) e\ {p+ po a\r{to)\0), 



(2.^ 



with (pq Pq a\J^{tQ)\0) the matrix element of the current creating the gluon from the vac- 
uum, with momentum {p^ = k~^ ,pq) and color a. Summation over repeated discrete indices 
is implied. Here, {k\ g^°'{tL,to) \po) = {kb\ g{tL,to) \poa) (we use indifferently both nota- 
tions) is the effective propagator that describes the (non relativistic) motion of the gluon 
in the transverse plane under the influence of the time dependent background field A~ . 
This propagation preserves k~^, as already mentioned, and Q depends on A;^, but we shall 
generally not indicate this dependence explicitly. Note also that the propagator Q do not 
carry any Lorentz indices. This is because, as shown in Appendix A, the propagation in 
the background held preserves the (transverse) polarization of the gluon. 

It is shown explicitly in Appendix A that Q is the propagator of a Schrodinger equation 
in two dimensions (the transverse plane) for a non relativistic particle of mass k^ moving 
in a time dependent potential A~(x^, x), with x'^ playing the role of the time. That is, it 
satisfies the equation 



W + 



2k+ 



{x\ g^^(x+, y+; k+) \y) = i6ab6{x+ - y+)5{x - y), 



(2.9) 



with iD = id + gA and d = d/dx^ . The solution can be written as a path integral 



{x\g{x+ ,y+;k+) \y) 



■ fc+ 



rx 
iy+ 



tj{x+,y^ 



(2.10) 



with r{y~^) = y, r{x~^] 
adjoint representation 



X, and f7 is a Wilson line evaluated along the path r{t) in the 



Uix+,y+-r) = Texp{ig / dt A" {t,r{t))T'' 



(2.11) 



Note that we work in a regime dominated by multiple scattering, so that the accumulated 
phase over a typical interval At is large, i.e., gA^At ~ 1. It follows that the exponential 
cannot be expanded in general (as done for instance in the opacity expansion [14, 16]). 

The probability to find the gluon with momentum k is obtained by taking the modulus 
squared of the amplitude in Eq. (2.8), summing over the final color indices, and then 
averaging over the color background field. Since the distribution of gauge fields is taken 
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Figure 1. Amplitude (upper part of the diagram) and complex conjugate amplitude (lower part) 
describing the propagation of a gluon in the medium from time to to time tL when it leaves the 
medium. The horizontal wavy lines represent the gluon propagator or its complex conjugate. Note 
that we use a bar to indicate color indices in the conjugate amplitude. At time t = t^, b = b. 
The (Gaussian) medium average generates instantaneous contractions that are depicted as vertical 
dotted lines, or as dotted circles. 

to be Gaussian, cf. Eq. (2.2), the medium averaging amounts to contracting pairs of fields 
from the Wilson hues in the two propagators representing the gluon in the direct and the 
complex conjugate amplitudes, respectively. These contractions can either connect the two 
gluon lines, or they can be 'tadpoles' with both endpoints on a same gluon line. It is 
convenient to represent the amplitude together with its complex conjugate on the same 
diagram, with time flowing from left to right. (This is somewhat analogous to the close 
time path technique, although we shall not explicitly use this particular formalism here.) 
The ensuing diagram, shown in Fig. 1, displays typical contractions. This representation 
makes it easier to discuss color and momentum flows. 

Let us consider first the color structure. At time tL, the system formed by the gluons 
in the amplitude and its complex conjugate is clearly in a color singlet state: both gluons 
carry the same color index b (i.e., 6 = 6 in Fig. 1), which is summed over since we do 
not observe the color of the final gluon. The field contractions involved in the medium 
averaging do not change the overall color state of the two gluon system, which remains 
therefore a color singlet at any intermediate time between to and tL- This implies that the 
medium averaged squared amplitude depicted in Fig. 1 contains a contribution of the form 
(as we shall see shortly k = k) 



■L,to)\po){Po\G^''\to,tL)\k)) = 6-^{k-k\S^^\tL,to)\po;Po) 



(2.12) 



where the sum over the repeated color index b is implied, and 



{k;k\S'''^'>{tL,to)\po;Po) 



1 



( 



Tj:{k\g{tL,tQ)\po){po\GHto.tL)\k)) , (2.13) 



where the trace Tr concerns color indices. The quantity S^'^\ referred to as a 2-point 
function, is the simplest of several n-point functions that we shall have to consider in this 



paper, and which are calculated exphcitly in Appendix B. Note that, in defining S^'^\ we 
choose to keep the time ordering given by the propagator in the amphtude, and to put all 
the momentum variables at the final time in the bra, and the ones at the initial time in 
the ket. The variables corresponding to the conjugate propagator are denoted by a bar, 
and are separated from the variables of the propagator by a semicolon. This notation will 
be used throughout (see also Appendix B). 

Similarly, since the correlator of the background field is invariant under translations 
in the transverse plane, the medium does not change the overall transverse momentum of 
the pair of gluons: whichever momentum is picked up from the medium by the gluon in 
the amplitude is compensated at the same time by a similar transfer to the gluon in the 
complex conjugate amplitude. Thus, if one chooses the final momentum to be the same 
in the amplitude and the complex conjugate amplitude, i.e., k = k {as we need to do if 
we consider the production of a gluon with a given transverse momentum) then, at each 
instant of time, the momenta remain equal in the amplitude and its complex conjugate. In 
coordinate space, as shown explicitly in Appendix B, this corresponds to the two gluons in 
the amplitude and its complex conjugate propagating with a fixed separation. The medium 
average of the Wilson lines in the pair of propagators is then of the form 

Cp{tL-to;r) = j^^--^(TTU{tL,to;x)U\tL,to-,x] 



exp { -g^NcitL - to) [7(0) - j{r)] } , 



exp 



-{tL-to)a{r) 



(2.14) 



where r = a; — ;r is the constant distance between the two gluons. This is recognized as 
the forward scattering amplitude for a color dipole which propagates through the medium. 
This is of course a fictitious dipole, which is formed with the gluon in the direct amplitude 
and that in the complex conjugate amplitude, which in the calculation of the probability 
group together in an overall color singlet state — the 'dipole'. The expression in the third 
line of the above equation, which features the dipole cross-section cr{r), will often be used 
to simplify writing in what follows. 

Altogether, the medium averaging of the squared amplitude will therefore generate a 
2-point function of the form (see Appendix B for details) 

{k;k\S^^\tL,to)\Po;Po) = (2^)^5(2) (po -po) V{k-po,tL-to), (2.15) 

where 

V{Ap,At) = j d^re-'^P-"- Cf\At;r), (2.16) 

can be interpreted as the probability that a gluon acquires a transverse momentum Ap 
while traversing the medium during a time At. This is best seen by using the harmonic 
approximation (2.6), where 

N^na{r)^]^qr'^ . (2.17) 
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Then we get 

(At;r) = exp(-^r2), (2.18) 

and 

P(Ap, At) = I dV e^r-Ar>-^^r^ ^ ^ (2.19) 

This formula confirms the physical interpretation of q: qAt is the transverse momentum 
squared acquired by a gluon after it has traveled through the medium for a duration At. In 
fact, in the harmonic approximation, 'P(Ap, At) solves the diffusion equation with initial 
condition V{Ap, At = 0) = (27r)^ 6^'^\Ap), with q playing the role of a diffusion coefficient. 
Note also that this probability is independent of the gluon longitudinal momentum . This 
is to be expected from the formal analogy of the present problem with the 2-dimensional 
diffusion of a non relativistic particle with mass k^. 

The exponential factor in Eq. (2.18) reflects a general feature of the interaction of color 
objects with the medium: the propagation of the color dipole (formed by the gluons in the 
amplitude and its complex conjugate) is not affected by the interaction as long as its size 
is small enough to be 'viewed' by the medium as a color singlet, but it is strongly damped 
as soon as the size of the dipole exceeds 2 /-/qAt. 

I I 
I I 

\ PO k \ 

to *L 

Figure 2. Diagrammatic illustration for the process in which a gluon of momentum {Pq,Po) is 
produced via a hard process (represented by the circled cross) and acquires transverse momentum 
k — po in its propagation though the medium (the thick wavy line). The + component of the 
momentum is conserved in this propagation, i.e., = Pq ■ 



We are now in a position to calculate the differential cross section to observe a gluon 
with transverse momentum k and arbitrary color and spin at time t^, after its production 
through some hard process at time tg. This is given by^ 

dao f d^po , , s dfJhard 



where dQk = (27r)~^d^fc dA;"''/2A;"'' is the invariant phase-space element, and do"hard/dr2pQ = 
I J"(po)P is the 'hard cross section' for producing the gluon (the sum over color a is under- 
stood - we also assume that all kinematical factors necessary to build a cross section are 



*The subscript on this cross-section refers to the fact that the process considered here involves no gluon 
branching. Similarly, the quantity ai in Eq. (2.21) below refers to a process which involves one branching. 
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included in the current J). In deriving this result, we have summed over the polarization 
vectors with the help of the completeness relation '^x^\i^)^*\ i^) ~ ■ This process, 
depicted in Fig. 2, has a probabilistic interpretation, with the wavy line representing the 
probability (2.19). 




^0 t 



Figure 3. Graphical illustration of the equation (2.21 ). The thick wavy lines represent the proba- 
bility V for transverse momentum broadening, the black dot is the splitting probability IC, and the 
circled cross is the cross section of the hard process producing a gluon of momentum po. 



It will be our main goal in this paper to show that the cross section for the process 
where one gluon splits into two gluons under the effect of medium interactions can be 
given an analogous probabilistic interpretation. Our subsequent calculations will lead to 
the following estimate for producing soft gluons {k^, k^ <^ ooc) 

X /^dt/ Vika-p,tL-t)Vih-q+p,tL-t) 

J to Jpo,q,p 

X /C(p - zq, z, p+) V{q-po,t- to) , (2.21) 

as 

and it is understood that z = k^ /pq. This result can be interpreted as a classical branch- 
ing process, illustrated in Fig. 3: after propagating from to to t, during which it acquires a 
transverse momentum q — po, the original gluon splits into gluons a and b with a probabil- 
ity ~ aslC{p — zq, z, g+) that depends upon the longitudinal momentum of the parent 
parton, the longitudinal momentum fraction z = jq^ carried by gluon a, and the trans- 
verse momentum difference p — zq. (The conservation of longitudinal momentum implies 
of course = = A;^ + k'^ with A;+ = p+ = zq^ .) After the splitting, the two gluons a 
and h propagate through the medium from t to t^,, and thus acquire some extra transverse 
momentum. 

Note that, in Eq. (2.21), the splitting occurs instantaneously at time t, that is, the 
effective splitting vertex /C(p — zq,z,q'^) is local in time. Moreover, the transverse mo- 
mentum is conserved at the splitting, meaning that one neglects the additional momentum 
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transferred from the medium to the gluons during the branching process. These are of 
course approximations, which are correct so long as the duration of the branching process, 
Tj^j,, is much shorter than the medium size^ L. Previously, in Eq. (1.2), we have estimated 
this time t^^^ for the case of asymmetric branchings, where one of the daughter gluons is 
much softer than the other one (say z <C 1). An estimate valid for arbitrary values of z 
can be obtained by considering AE, the change in the light -cone energy (the minus com- 
ponent of the 4-momentum, or equivalently the change in the energy in the equivalent non 
relativistic 2-dimensional problem), at the splitting vertex: 

AE=^ + A^^lP)^ - ^ = JPz^ . (2.22) 
2p+ 2{q+-p+) 2q+ 2z{l - z)q+ ^ ^ 

The inverse of this energy is the formation time, which gives access to the duration of 
the branching process. To get this, all what one needs to do is, as in Eq. (1.2), relate 
p — zq to the transverse momentum acquired by multiple scattering during t^^^. That is, 
(p ~ zq)"^ ~ ^■^br' which together with the above expression for AE' = 1/tj^^ provides an 
estimate for the corresponding branching time: 



z(l-z)p+ 

-. = ^^-y^, (2.23) 

which generalizes Eq. (1.2) by the replacement of w = zp'^ by z{l — z)p'^ (note that the 
latter quantity can be interpreted as the reduced mass for the effective non-relativistic 
motion of the two produced gluons in the transverse plane). This makes it clear that 
Eq. (1.2) holds only for asymmetric splitting with z <C 1. A slightly more precise estimate 
for r^^ will be given in Sect. 5 (Eq. (5.5)). 



3 General structure of the in-medium gluon branching 

In this section we analyze the main structure of the calculation of the gluon branching, and 
show how it can be conveniently divided into three main stages: (i) production of the gluon 
at time to via some local (hard) process, followed by the propagation of the gluon in the 
medium; (ii) the process where the gluon splits into two gluons; (iii) finally the propagation 
of the two new gluons until the end of the medium at coordinate = tL = \f2L. 

3.1 The amplitude 

The amplitude for the branching process is illustrated in Fig. 4: a gluon is created at time 
to with color co and momentum (p^,Po) out of the current J'^^^p^ ^Vo) = (Po co|'/(to)|0); it 
then propagates in the medium where it acquires a transverse momentum Qi—po, changes 
color form cq to ci, and eventually splits into two gluons at time ti. The newly produced 
gluons then propagates through the medium until the time t^, where they are observed 
with momenta ka,kb, colors a,b and polarization Xa, ^b- Using the elements described 

^In the presence of multiple emissions, the correct reference scale for T^^^ is not L, but the typical time 
interval between two successive emissions. We shall return to this issue in the concluding section. 
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tn ti tr 



Figure 4. Graphical illustration of the amplitude for one gluon splitting. The initial gluon produced 
at time to with momentum po s^nd color cq splits at time ti into two gluons with momenta pi and 
P'l (= Qi ^ Pi) respectively. Note that the + component of the momenta is conserved both in the 
propagation and in the splitting, so that k'^ = zpq , fcj^ = (1 — z)p^ . 

in Appendix A, together with compact, matrix, notations similar to those already in the 
previous section (cf. Eq. (2.8)) one easily gets the following expression for this amplitude: 

Mt,M ka, K h) = r 

X / dti{kaa;kbb\g{tL,ti)g{tL,ti)\piai;p[bi) 

J to 

X (pi ai;p[ bi\r^^ \qi ci)(qi ci| ^(ti, to)|Po cq) J'''°(Po 'Po) , 

(3.1) 

with = fc^/(2A;+), k^ = k'^/{2k^) and z = /pq . A summation over repeated discrete 
indices is implied. We have defined 

{ka a; kb b\ G{tLM)Q{tLM) \Pi ai; p'l bi) = (fcal^'^"^ (^L, ti) \Pi){kh\ G^^' (tL, ti) Ip'i) , 
(pi a,-p[ b,\r^' ki ci) = i27rf5ipi+p[ - Qi) 2g f^'^'^^r^^pi - zq^z) , (3.2) 

and we denote indifferently the matrix elements of the single propagator by {kalG""""^ (^L) ^i) |Pi) 
or by {kaa\G{tL,ti) \pi oi), and similarly (pi ai;p[ bi\ P^' \qi a) = T^ffc^ci l^i)- 

To calculate the cross section, we have to multiply the amplitude by its complex 
conjugate amplitude. In order to analyze the effect of the medium average, it is convenient 
to draw the amplitude and the conjugate amplitude in the same diagram, one just below 
the other as shown in Fig. 5 (and in Fig. 1 of the previous section). The upper part of 
the diagram corresponds to the amplitude and the lower part corresponds to the conjugate 
amplitude, with time flowing from left to right identically in the amplitude and its complex 
conjugate. Gluons which are drawn parallel to each other represent the same physical 
gluon (in the amplitude and conjugate amplitude, respectively), so they have the same 
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Figure 5. The three regions that are involved in the description of a gluon sphtting in a medium. 
Time flows from left to right and the amplitude (upper part) is drawn together with the complex 
conjugate amplitude (lower part). The + momentum is conserved in the splitting. The various 
indices at each ends of the gluon lines (wavy lines) are color indices. The times separating the 
various regions are light-cone times: ti and t2 are the light-cone times where the gluon splitting 
occurs in the amplitude and the complex conjugate amplitude, respectively. The indices a, 6, c are 
color indices, with each gluon carrying a given set of color indices that involve the same letter: thus 
'gluon a' can be seen in the various colors a, ai, 02, etc. 



momenta and same color at t = t^. For definiteness, we have chosen the splitting to 
occur first in the amplitude and then in the conjugate amplitude (i.e., ti < t2), the inverse 
ordering can be accounted for by taking 2 times the real part when calculating the cross 
section (see Eq. (3.13) below). One may then identify three distinct time regions that 
correspond respectively to the propagation between to and ti (region I), between ti and 
t2 (region II), and between t2 and t^ (region III). To facilitate the grouping of terms in 
the different regions, we make repeated use of the identity (A. 20), writing for instance (in 
matrix notations) G{tL,ti) = G{tL,t2)G{t2,ti), or ^"''(to,t2) = ^"'"(to, ti)^"''(ti, t2)- Next, 
one introduces closure relations at appropriate times, over complete set of states in the 
transverse plane and in color space. The matrix element in the amplitude can then be 
written as (see Figs. 5 and 7 for pictorial representations) 

(fc,;fc6|g'^'^Hti,t2)a"^'^Hi2,tl)g^^^(ti,t2)g^^^K*2,tl)rif,^ei S''''ih,to)\po) 

= [ (fca|g'^"^(tL,t2)|q2)(q2|a'^^'^Hi2,tl)|pi)(fc5|a''^(tL,t2)|g^)(q^|a'^'Hi2,tl)|p'l) 
■J <l2,Pl,p'i,Q2 

x{pr,p\\T%^Jq,){q,\g^^^^{h,to)\po), (3.3) 
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and that in the complex conjugate amphtude as 

''Q1,Q2,P2,P'2 

x{p2\g^~^'%t2,tL)\ka){p'2\g^'^\t2,tL)\h). (3.4) 

One can now multiply the matrix elements of the amplitude and its complex conjugate, 
and regroup terms in the various regions. One writes the result in the following way 

ika\g''''HtLM\Q2)ikb\g''HtLM\Q2)iP2\g^~'''^^ 

x{q2\g''^''Ht2,ti)\pi){q'2\g'''Ht2,ti)\p\){q^^^^ 
>^{qi\g'''"{ti,to)\Po){Po\g^'''Hto,h)\q,), 



(3.5) 



where the three lines in the expression above correspond to the three regions III, II, I, 
respectively, from top to bottom. At this point, we are ready to discuss the medium 
average and identify the simple, factorized, color structure of the three regions. 

3.2 Color structure 

As already noted in the previous section, the contractions of background fields that are 
involved in the medium average preserve the overall color of the system, and since all such 
contractions are instantaneous, the overall color state of the system comprising all gluons 
in the amplitude and its complex conjugate is conserved at all times. Since, as we shall 
verify shortly, the overall color state is a color singlet at t = ti, the gluons traversing the 
medium are in an overall color singlet at any given time. This property allows us to perform 
the medium averages explicitly by reducing considerably the number of color states to be 
considered along the calculation: all one needs to do is to find the appropriate singlet state 
in which the system is at the separation points ti and t2- The resulting color structure is 
illustrated in Fig. 6. 

That the system of gluons is in a color singlet state at t = tj;, is easy to see. Indeed the 
colors of the gluons are identical in the amplitude and in the complex conjugate amplitude, 
and since we do not observe the colors of the produced gluons, these colors are summed 
over, projecting the two pairs of gluons a and b onto on color singlets. Since the system 
remains a color singlet at all time, it is a singlet at time to, corresponding to the beginning 
of region I, with which we start our detailed analysis. 

Consider then region I, that is, the third line of Eq. (3.5). As we have just argued, 
the gluon system is in a color singlet state. Hence at time to, the gluon carries the same 
color index in the direct and complex conjugate amplitude (cq = cq), and then the medium 
averaging yields 



5^°^"°((qi|a^^'=°(ti,to)|Po)(j5o|at'°'Hio,ti)|gi)) =5'^'nqi;gi|5^'Hii,io)|Po;?5o), (3.6) 
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(I) 




{II) 




{III) 




Figure 6. The color connections that are responsible for making the system of gluons a color singlet 
in each of the three regions. The dotted vertices represents / symbols, /"^^, while the 'vertical' 
connections among them represents Kronecker delta's, dab- 



where the quantity S{'^\ti,to) has already been introduced in Eq. (2.13) and evaluated in 
Eqs. (2.15)-(2.16). Note that in Eq. (3.6) a sum over cq is performed. However, this sum 
is not free since the current also carries the same color index. However, since each term of 
the sum is eventually independent of cq, we can correct the result by simply dividing by 
— 1. In summary, the contribution of region I after medium average is 

(qi;qi|5(2)(ti,to)|Po;Po). (3.7) 



Now we turn our attention to region H, that is the second line of Eq. (3.5). We can use 
the Kronecker delta d^^'^^ from Eq. (3.7) and recall that the color structure of the vertex 
at ti is of the form to realize that the relevant correlator for this region is 

r'^^^((q2|a"^"Hi2,ti)|pi)(q^|a'^'Ht2,ti)|p'i)(qi|at^^^"^(ti,t2)|g2)>, (3.8) 

from where it is clear that we have a color singlet state at ti by contracting the three initial 
color indices with an / symbol. This color structure is unchanged by the medium average 
in this region given that there is no other way to combine these three gluons into a singlet 
state. It follows therefore that the color indices at t2 must also be contracted with an / 
symbol. More specifically, this means that Eq. (3.8) takes the form 

r^'i'=^((q2|a"^"Ht2,ti)|pi)(q^|a''^''Hi2,ti)|p'i)(qi|at^^^"Hii,i2)|q2); 
^ r''''Hq2q'2;q2\S'^'Ht2,ti)\pip[;qi), (3.9) 
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with 



{q2q2;q2\s^'^\t2,h)\pip[;qi) 



Nc{N^ - 1 



(3.10) 



Due to the simple color structure of this object, one can calculate explicitly the medium 
average of the Wilson lines inside the propagators, by resumming the contractions of each 
pairs of gluons independently. One gets 



= exp|-^^^ dt [a{rb - Vc) + a{rb - ra) + a{ra - Vc)]^ . (3.11) 

The color factor in front of the leading order dipole cross section can be calculated from the 
appropriate contraction of the color factors in the adjoint representation f"-bcjcdejefa _ 
— The calculation of the path integral needed to complete the evaluation of S^'^^ is 

done in Appendix B.2. 

Turning our attention to region III, the first line of Eq. (3.5), we see that the four 
gluons form a singlet state through the combination of two /'s, one coming from Eq. (3.9) 
and the other from the vertex at t2 (see Eq. (3.5)). As opposed to the previous cases, there 
are several ways of combining four gluons to form a singlet state. The final color state is 
given by the identification of two pairs of gluons as the conjugate of each other as explicitly 
shown by setting the color indices equal for the respective propagators at ti (as already 
discussed, a = a, b = b). The relevant correlator for region III is therefore 

r^>'^^^f^'^^^^{ka\g''-^tL,t2)\q2){h\g''H^ 

= NciN^ - l){kakb;kakb\S^^\tL,t2)\q2q2;P2P2)- (3-12) 

The fact that there is no unique way to combine the four gluons into color singlets, in 
contrast to what happened for the 2 and 3-point functions, prevents us to give an explicit 
expression for the average of the Wilson lines in the propagators analogous to Eq. (3.11) 
for the 3-point function. A more detailed analysis is required, which will be carried out in 
section 4. 

3.3 Momentum structure 

The fact that the correlator of the random background field depends only on the difference 
of the transverse coordinates leads to total transverse momentum conservation. It follows 
that at each instant of time, the sum of momenta in the amplitude equals that of momenta 
in the complex conjugate amplitude, since this is the case at the final time tL. As explicitly 
shown in Appendix B, each of the n-point function contains a corresponding (5-function 
expressing this property. There are also additional (^-functions hidden in the vertices F 
which express the conservation of transverse momentum in the local splitting. Using all 
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Figure 7. Flow of momenta in the amplitude for one gluon splitting and its complex conjugate. At 
any given time, the sum of momenta in the amplitude equals the sum of momenta in the complex 
conjugate amplitude. The + component of the momentum is conserved, k'^ + k'^ = Pq , and is not 
indicated. 



these (5-functions, one can then reduce considerably the number of momentum integrations 
when building the cross section for the branching process. The resulting flow of momenta is 
displayed in Fig. 7, which also lists the independent momenta that we use. We denote with 
a tilde the n-point functions from which the delta function for momentum conservation 
has been factored out. 

There are further simplifications that reduce the number of independent momentum 
variables. For instance, the 2-point function depends only on the difference of two momenta. 
Similarly for the 3-point function. Also the vertex depends only on p — zg, as visible on 
Eq. (A. 23). These extra simplifications will be exploited when needed. 

We are now in a position where we can write the cross section for the two gluon 
production process. After performing the medium average and summing over polarizations, 
one obtains 



J to J to J Po,Piqiq2P2q2 

92,92 -q2]P2,q2 -P2){q2,q2 - q2;q2\S^^\t2,ti)\pi,qi -pwqi) 

x(qi;qi|S(2)(ti,to)|po;Po)J^''=°(p^,Po)J*^'^»(p^,Po). (3.13) 



X 



In this equation, we used the compact notation dfi^ = (27r) ^d^A; dA;^/2A;+ for the element 
in phase-space and in the r.h.s. it is understood that pQ=k^ + and z = /pq . 
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3.4 Qualitative comments 

Before getting into the details of the evaluation of the various n-point functions that enter 
the expression above, we find it useful to make a few qualitative observations that will 
lead us to better appreciate the physical content of the calculations to be performed in 
subsequent sections. 

As we have seen, the n-point functions that appear in the various regions contain a 
medium average of a product of Wilson lines in an overall color singlet state, inside a 
path integration. The way to proceed in the explicit average evaluations is to calculate the 
average of a particular combination of Wilson lines with fixed paths in transverse coordinate 
space, then perform the integration over these paths. In general, the medium average of the 
Wilson lines decays exponentially with the distance between particles constituting a color 
neutral state, typically as Aa;^ > 2/(gAt) (see e.g. Eq. (2.18) for the 2-point function). 
That is, one gets a significant suppression factor whenever the particles are farther apart 
than the inverse of the typical transverse momentum acquired by multiple collisions. On 
the other hand, in the path integral, the gluons trajectories diffuse in transverse coordinate 
space, which typically increases the average distances between the gluons so that Aa;^ ~ 
At/cj. Combining these two estimates, one gets that in order avoid a large exponential 
suppression factor the longitudinal extent of the medium average in consideration should 
be hmited to At < ^J^hJfq = r^^ < L. 

The argument above does not always hold since it assumes that the diffusion in trans- 
verse space is independent for all the gluons, so that the mean distance between them is 
monotonically increasing. However, when two gluons are conjugate to each other their 
diffusions in transverse coordinate space are correlated, indeed their transverse separation 
stays constant (see Appendix B). This implies that the longitudinal extent of region I is 
not constrained by the argument above. On the other hand, the three gluons in region II 
propagate independently. This, combined with the fact that there is only one way to form 
a color neutral state with three gluons, implies that this region II must be short-lived, i.e., 
*2 — ^1 ^ ''"br! ™ order to avoid suppression factors. 

In region III we run again into the case where gluons can be paired with their respec- 
tive conjugates and therefore do not propagate all independently. Moreover, by pairing the 
gluons this way, two color neutral subsystems are formed, which can propagate indepen- 
dently and be arbitrarily far away from each other without introducing suppression factors. 
Nevertheless, the general argument concerning the suppression of correlations by medium 
rescattering still holds within a narrow region bordering region II, where all four gluons are 
still correlated. We shall find that the longitudinal extent of that intermediate region is of 
order t^^^ for the same basic reasons as for region II. Accordingly, we shall conclude that 
the four-point correlator obtained for the medium average in region III is factorizable into 
two two-point correlators except for a parametrically small (~ t^JL) region right after the 
second splitting. 
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4 Factorization of two-gluon propagation 



We turn now to the detailed analysis of region III. As already anticipated in the previous 
section, the main goal here is to show that the propagation of the two resulting gluons can 
be considered as independent and correlations among them are a subleading effect. In that 
case, all the effects of the splitting can be included inside the splitting kernel of region II, 
which will be calculated in detail in the next section. 

We need to calculate explicitly the 4-point function S^^\tL,t2)- Since the definition 
of the scalar propagators entering the average in (3.12) take a simpler form in coordinate 
space, it is convenient to calculate first the medium average of the four propagators in 
coordinate space with arbitrary endpoints, perform the Fourier transform, and then impose 
the necessary constraints in the momentum variables. 

The coordinate space 4-point function under consideration takes the explicit form 

{yayb;yayb\s''^\tL,t2)\ 

= j VraVnVfaVn expS^^j^'^dt {ktrl + k+rl - k+fl - k+fl) | 

^ja,b,c,j-a,h-c, (^lJ^^^^tL,t2; TaMb, {tl, n)Ul,{t2 , t f a)Ul,{t2 , t n)) , (4.1) 

where ra^b{t2) = Xa^bifafiiti) = ya,b and similarly for the bar coordinates. 

We start by first calculating the medium average of the Wilson lines before attempting 
to perform the path integrations. Since for the moment we are only concerned with what 
occurs in region III and all the propagators have the same longitudinal extent, longitudinal 
coordinates will be dropped from the Wilson lines, and only the transverse trajectories will 
be kept explicit. 

The calculation of the correlator in the last line of Eq. (4.1) turns out to be extremely 
complex even for our simple, Gaussian, model for the random background field. So, before 
we start the calculation, it is useful to make some remarks about its general structure and 
the strategy that will be employed for the rest of the section. Since our goal in this section 
is to show that, to the level of desired accuracy, one can factorize this correlator into the 
propagation of two independent gluons, it will be convenient to split the evaluation into a 
'factorizable' piece and a 'non-factorizable' piece where, 

S'-'HtL,t2) = Si2{tL,t2) + SSL(tL,t2), (4.2) 

with 

(Z/a2/fe; yayb\Sf2iiL, t2)\XaXh; Xa, X^) 

= {ya]ya\S''^\tL,t2)\Xa]Xa){yb-,yb\S^^HtL,t2)\Xb-,Xf,), (4.3) 

or, in terms of the momentum representation entering the expression for the cross section, 

{kah; kah\S^2(.^L,t2)\q2, q2 - q2;p2,q2 - P2) 

= {27rf6'^^\p2-q2)V{ka-q2,tL-t2)V{kb-q2+q2,tL-t2). (4.4) 



- 21 - 




^2 tr to tr 



Figure 8. Graphical representation of the correlator (4.6) and its large-A'c version (4.7). The 
horizontal lines in between the vertical dashed lines represent the Wilson lines either in the adjoint 
(gluon lines) or fundamental (quark lines) representation. Lines in the right and left ends of the 
diagrams indicate color connections only, with gluon lines identiiying adjoint color indices, quark 
lines identifying fundamental color indices, and three gluon vertices representing / symbols. 



The analysis below is therefore aimed to show explicitly that the non-factorizable pieces 
are suppressed with respect to the factorizable pieces and therefore it is only the expression 
in (4.4) which enters the result for the cross section of the two-gluon production process. 

The complete evaluation of the correlator in (4.1) is complicated because of the many 
singlet states that can be constructed with four gluons, but by transforming this expression 
into one with only Wilson lines in the fundamental representation it is possible to gain some 
insight about the transitions between the different color states. In order to do that, one 
makes use of the following identity relating Wilson lines in the adjoint representation with 
Wilson lines in the fundamental representation. 



Uab{r) = 2TV C/t(r)f'^C/(r)t' 



(4.5) 



After repeated use of this identity, one can perform the color algebra and eliminate all the 
explicit color factors (see details in Appendix C), arriving at 

r^'^'V'^'^'^(?:>aa,(ra)[/6fe,(r,)[7i^,(f,)c7i^(n)) 



-Tr 



U{ra)UHfa)\ Tr \u{n)U\ 



Tr 



u\ra)U{n)u\n)u{fa] 



u{ra)u\fa)u{n)u\n)u{fa)u\ra)u{n)u\n) 



+ h.c. 



(4.6) 



This may not look like a real improvement: the number of Wilson lines has doubled 
and now we have several nontrivial terms to deal with! The advantage of this approach 
is that it greatly simplifies when one considers the large- A'^c limit, a common approach 
used to simplify the calculations while capturing the essential physical behavior. In the 
fundamental representation, each color trace is proportional to Ac, from where it becomes 
clear that the first term will dominate over the second in the large-Ac limit. Still in 
this limit, one can factorize the average of a product of traces as the product of the 
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independent averages of each individual trace. Indeed, correlations between fields entering 
different traces are suppressed by inverse powers of Nc- By also using the fact that all the 
considered correlators are real (at least, in our Gaussian model for the medium averages), 
we finally get the following expression for the dominant contribution at large Nc : 



Tr 



Tr 



n 



Tr 



U^ira)Uin)U\n)U{fa) ). (4.7) 



This relation can be easily seen graphically by replacing the gluons with quark-antiquark 
pairs as shown in Fig. 8, where only the planar part was kept. Each of the loops corresponds 
to one of the traces in (4.7) and connections among them are suppressed. 

The three factors in (4.7) can then be evaluated independently. Two of them are quark 
dipoles which take the by now standard form (compare with Eq. (2.14)) 



1 



Tr 



U{r)U^f) 



exp 



ftL 

/ dta{r — r) 

Jt2 



(4.8) 



The third factor in (4.7) has the form of a quadrupole amplitude. This kind of cor- 
relator has been studied in the literature, always for Wilson lines with fixed transverse 
coordinates while here we need to be able to evaluate it for arbitrary trajectories. A 
straightforward generalization, as it was possible in the case of the dipole amplitude, is 
so far not available. On the other hand, the use of the large-A^c limit helps greatly to 
further simplify the calculation and allows us to gain some insight into which terms are the 
important ones for the situation at hand. 

There are only two independent singlet states that can be formed with two quark- 
antiquark pairs corresponding to the two possible ways of pairing quarks with antiquarks. 
Considering the way the color indices are contracted for the quadrupole amplitude, one can 
consider this quadrupole amplitude as the probability of switching from one singlet state 
to the other one, after undergoing multiple scatterings with the medium. For the case at 
hand, it is clear that the natural causal way to consider the process is such that the two 
singlet states are made with the pairs of lines o, b and respectively a, h in the initial state, 
and with the pairs a, a and 6, b in the final state. 

The two aforementioned singlet states are not orthogonal and therefore one can isolate 
a contribution where the system was initially in the final color state and all the interactions 
with the medium preserve this color state. In that case, one can consider the four-particle 
system as two separate quark dipoles, which do not see each other, yielding a product of 

(2) 

two already familiar dipole amplitudes Cq . This type of contribution has the exact same 
structure as what was referred at the beginning of the section as the factorizable piece and 
the calculation below will show explicitly how one arrives at Eq. (4.3). 

The non-factorizable piece consists of all the possible cases not accounted for in the 
factorizable piece, which are those where transitions between the color states are allowed 
during the interaction. In general one can have an arbitrary number of such transitions, but 
in the large- A^c limit those are suppressed by inverse powers of Nc- In that case, one must 
take into account only contributions where there is only one transition between the two 
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Figure 9. Sample of the kind of interactions allowed at each end of the quadrupole amplitude. 
The dashed line in the middle represents the transition point in between the two regions where a 
different kind of contraction is allowed. 

possible singlet states. All the diagrams which are not suppressed by inverse powers of Nc 
have the general form depicted in Fig. 9: on the left side of the diagram only contractions 
between lines a and b, or between lines a and b are allowed, while on the right side of the 
diagram only contractions between the lines a and a, or between lines b and b are allowed. 
These two regions are separated by just one contraction performing the transition, which 
can connect either lines a and b, or lines a and b. It is easy to show that those are the only 
planar diagrams present in this calculation. 

In order to find a closed expression for the sum of all these diagrams one has to choose 
an orientation to perform the resummation. For our purpose, it is convenient to read the 
diagrams from right to left. In that case, one starts with the dipole formed by lines a, a 
and b,b up to ^3 at which point the transition contraction occurs. Specifically, at ts, a 
transition occurs to the singlet state where the two dipoles are formed by the pairs a, b and 
a, b. All the contractions on both sides of the interactions can be resummed into dipole 
amplitudes over the corresponding longitudinal extent. This procedure assumes that there 
is a transition, and therefore ignores all the diagrams with no transition, i.e., the diagrams 
where all the contractions respect the singlet structure in the right side of the diagram. In 
order to properly account for those diagrams one must include the factorizable contribution, 
given by the product of two dipole amplitudes for the respective dipoles in the right hand 
side of the diagram. The total result is then the following: 
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Figure 10. Contractions contributing to the transition amplitude between the two singlet states. 



QitL,t2;ra,n,ra,rb) 



Tr 



u\ra)Uin)uHn)uifa) 



+ 



tL 



xCf )(t3, t2; r-a, n)Cf'\t^M\ ra, n), (4.9) 



t2 



where ^(ta) denotes the transition between the two states. It receives contributions from 
the diagrams in Fig. 10 and is explicitly given by 



T{t3) = Cpn [a{ra - n) + a{ra - n) - a{ra - n) - a{ra - n)] 



*3' 



(4.10) 



where all coordinates are evaluated at ^3. The separation between the factorizable piece 
and the non- factorizable one is manifest in formula (4.9). It is important to note that the 
same analysis can be performed by reading the diagrams from left to right. One obtains 
a formula which is the same as above except that the labels b and a are interchanged. In 
that case the transition amplitude includes a, a and 6, h contractions instead of the first two 
diagrams of Fig. 10. It is easy to show that the two expressions are equivalent by noticing 
that if one substracts the two different ways of writing the non-factorizable term one can 
perform the is integral explicitly and show that the resulting contribution vanishes. 

Now we are in a position to insert the medium average of the Wilson lines into the 
4-point function (4.7). Consider first the term corresponding to the factorizable piece of 
the quadrupole, one gets 



VaVb] VaVb •^a-^b, "^ai "^b 

VraVnVfaVfb exp 



t2 



Cf\tLM]ra,fa)Cf\tL,t2; 



2 



(4.11) 



In the large- A'^c hmit the square of a quark dipole amplitude is equivalent to a gluon dipole 
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amplitude, therefore this factorizable piece can be rewritten as 



{yaVb; VaVblSf^itL, t2)\XaXb; Xa, Xb) 



= j VvaVfaexpl^^ J^"^ dt {rl-f-l)^CP{tL,t2;ra,fa) 

X j VnVfb exp I ^ ^ ^ dt {rl - fl) | Cf\tLM] n, fb), 

= {ya;ya\S^'^HtL,t2)\xa]Xa){yb;yb\s'-^HtL,t2)\xb-,Xb), (4.12) 

where we have identified the path-integral representations of the 2-point function S^^^ (ti, 12) 
according to Eq. (B.4). We thus find the factorized structure anticipated in Eq. (4.3), at 
least within the large- A^c approximation. 

Now let us focus on the non-factorizable piece. The last line on (4.9) has the same 
structure as the factorizable piece except for the fact that its starting longitudinal coordi- 
nate is at ^3. One can still recombine this piece with the remaining dipole amplitudes in 
(4.7) since those amplitudes can be split in two in the following way: 

C(2)(t^,i2;r,r) = C^^\tL,U;r,f)C^^\h,t2;r,f), (4.13) 

yielding to two independent gluon dipoles for the region with longitudinal coordinate 
greater than ^3. 

{yayb\ yayb\S^l^{tL, t2)\XaXb; Xa, Xb) 

= J VraVnVfaVfb expl^^J^^'dt [ktrl + k^rl - k^^l - k+f-^) | 



L 



dts Cjj^HtL,t3;ra,fa)C^'^HtL,t3;rb,fb) T{t3) 

xCf ) (t3, t2;ra, n)Ci'^^ (ts, ^2; r^, n)Cf ) (t3, t2;ra, rb)Cf^ {h, t2;ra, n)- (4.14) 

Moreover, in the same way that the scalar propagators follow the convolution relation 
(A. 20), one can also split the longitudinal extent of the path integrals at and explicitly 
perform the integration in the region where the system factorizes into two independent 
gluons. One gets, 

{yayb;yayb\sl2,c(^L,t2)\xaXb;xa,Xb) 



2 



r dt3 I _ _{ya;ya\S'^^\tL,t3)\za;za){yb;yb\s'^^\tL,t3)\zb;zb)T{Z) 
X ^ ' VraVnVfaVfb expS^'^l^'dt {k+fl + k+f^ - k+fl - k+fl) | 



'i2 I ^ Jt2 

xCf )(t3, t2; r-a, rb)Cf\t3M] Ta, n)Cf )(*3, ^a, n)Cf )(*3, t2] fa, n), (4.15) 

with rafi{t3) = Zafi and similarly for the bar coordinates. 

The transition amplitude T depends explicitly on the new intermediate coordinates. 
One can easily see that in the harmonic approximation it takes the simple form 

T{Z) = -Uza-Za)-{Zb-Zb). (4.16) 
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When switching to momentum space, these coordinate differences can be expressed in terms 
of derivatives of the two-point functions going from ^3 to tL- 

After performing the Fourier transform, the non-factorizable piece can be written as 



{kah]kakb\Sl2,citL,t2)\q2,q2 - Q2;P2,q2 -P2) 

/ dts / I VP(fca - qsa, tL " ^s) " VP(fct, - qsb, tL " h)i{Q3a, Qsb, Q2, Q2, P2), 

(4.17) 



where V denotes the gradient with respect to the transverse momentum variables, and X is 
the Fourier transform of the path integral in (4.15) after removing the overall momentum 
conserving (5-function, 

(2^)2<5(2)(g^ _ + q2)i(Q3a, Qsfe, ^2, Q2, P2) 

i[qSa-(Za-Za) + q3biZb-Zi,)-q2-Xa-q'2-Xt+P2-Xa + {q2-P2)-Xi,] 



e 

{x,z} 



X / VraVnVfaVn exp <j - / dt [k+fi + k+rt - k^K - k, 



^^a ' '"^^^ '"^^^ 

t2 Jt2 

xCf )(t3, ^2; r^, n)cf )(t3, t2;fa, n)cf\h, t2;n, n)cf\h,t2;ra: n). (4.18) 



H)] 



This object is explicitly calculated in Appendix B.3. The main result of that calculation 
is that the result of the path integration is not exponentially suppressed only if its lon- 
gitudinal extend is small, of order of the formation time. In terms of our result for the 
non-factorizable piece of the 4-point function it means that we can restrict the integration 
over t3 to only the interval between t2 and t2 + tf, leading to the the fact that the region 
after is of the order of the length of the medium while the region between t2 and is of 
the order of the formation time. 

Using the explicit expression for the "P's in Eq. (2.19), one can easily see that the 
derivatives in (4.17) give a factor of 

(fcg -qsa)- {kb - qsb) J_ . 
(PitL-hY ^ qL- ^ ■ 

This additional factor of L in the denominator can not be compensated in any way since 
the integration in t3 was shown to have support over a small region and therefore does not 
depend on L. Therefore, the non-factorizable piece of 5"^^^ is parametrically smaller than 
the factorizable piece by a factor of t^J L and can be safely discarded for the kinematical 
regime under study. 



5 The splitting kernel and the gluon-splitting cross section 

In this section we shall complete the calculation of the cross-section for medium-induced 
gluon branching, establish the formula (2.21), and give an explicit expression for the split- 
ting kernel /C. 
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We return to Eq. (3.13), and make a first simplification that exploits the main result of 
section 4, namely the factorization of the 4-point function. Thus, we replace in Eq. (3.13), 



S^'^\tL,t2) by s\2iiL,t2) given by Eq. (4.4). We obtain then 



23fie 



tL 



dti 



to 



ti 



dt2 



to 



xr^'{pi-zqi,z)r'^'{p2-zq2,z) 

xV{ka - 92, - t2)V{h - 92 + 92, - ^2) 

x(p2,q2 - 92; 9215*^^^*2, ii)|Pi, 91 -Pi;9i) 
xP(9i -Po,ti -*o)>'"(p^,Po)J"''°(p^,Po 



(2^)V2)(p2-92) 



Po,pi,P2,qi,q2q2 



(5.1) 



where = A;,^ + /c^^, 2 = k'^/pQ. (In order to follow the flow of momenta it may be 
useful to refer to Fig. 7.) At this point we recognizes in Eq. (5.1) the three 'classical 
propagators' V expressing the transverse momentum broadening of the gluons before and 
after the branching. These are almost the same as in Eq. (2.21), except that the momenta 
are those appropriate to regions I and III, while in Eq. (2.21) the finite extent of region II 
is neglected. We shall return to this question shortly, and identify now the splitting kernel 
/C by combining the 3-point function S^^\t2,ti) with the vertex factors that are explicit in 
Eq. (5.1). 

Consider first the two vertex functions. They combine to yield 



r^\Pi,z)r^\P2,z) 



1 



2 + (1-Z)2 



Pi- P2S'' + 2PIP^, 



(5.2) 



where we have set Pi = pi — zqi, P2 = P2 — zq2- For simplicity, we consider here only 
inclusive cross sections that are averaged over azimuthal angles. Under this assumption 
one can then replace 2 PI by (Pi • P2) (5**, and get 

4 



where 



NcV^\Pi,z)r'j\P2,. 



Pgg{z) = N, 



zl 



Pi • P2) 5' 



(5.3) 



+ z{l-z) 



is the leading-order Altarelli-Parisi splitting function [37]. 

Next we consider the 3-point function S^^\t2i ti), which is calculated in Appendix B.2. 
In the harmonic approximation it reads (cf. Eq. (B.28)) 



{p2,q2-P2;q2\S^^\t2M)\ 
87r[l + z2 + (l 



Pi, 91 -Pi;9i) = 

2[l + z^ + {l-zmqi-q2f 



3qAt 

27r(l + i) 
kl sinh(OAt) 



exp 



exp 



n + i) 



3qAt 

{P1+P2? 
Akl coth(17At/2) ' 



(1 + 



(A - P2? 1 
4/c2 tanh(0At/2) [ ' 



(5.4) 
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with At = t2 — t2, and 



1 + i z(l — z)p'^ 2 l~t I 

^ = ^— > -^b. = y ^ > = ^eff r^, = V z(l - ^)Po 9eff , (5.5) 

br V yen 

where ^eff denotes an average, ^-dependent, version of the jet quenching parameter: 

gcff = \q [z^ + (1 - zf + 1] . (5.6) 

In writing the expression above for the 3-point function, we have used the delta-function 
5^'^\p2 — Q2) in Eq. (5.1) in order to replace q2 by p2- This allowed us in particular to 
express S^^\t2,ti) in terms of -P2. This expression (5.4) of S^^\t2,ti) makes it clear that 
the extent of region II is limited: indeed, since l/sinh(|r2|At) oc exp{-At/\/2rbJ for large 
At, the time separation At = t2 — ti needs to be kept smaller that t^^ in order to avoid 
exponential suppression. Thus, as anticipated via simple considerations in Sect. 2, and in 
the qualitative discussion of Sect. 3.4, one sees that the time scale of the branching process 
is indeed determined by t^^. Recall that this is also (approximately) the same time scale 
that controls the exponential damping of the non factorizable part of the 4-point function 
(see Sect. 4 and appendix B.2, Eq. (B.51)). 

In order to proceed further, we shall exploit the fact that region II has limited extent 
and neglect in the factors V the fraction of momenta that can be attributed to momentum 
broadening within region II. Since this region is of extent r^^., this fraction is typically of 
order t^^JL as compared to the momentum acquired through propagation before and after 
the splitting (this is in fact explicit for the dependence of S^^-* on qi — q2 in Eq. (5.4), where 
it can be seen that {qi — 92)^ is at most of order qT^^)- Thus we shall replace in Eq. (5.1) 
V{kb-q2 + q2,tL-t2) by V{kb-qi + q2,tL-t2) (ignoring the small difference gi - 92 ) , and 
use as independent variables qi, Pi, P2, qi — 92 in place of Qi , Pi , Q2 j Q2 • Since the factors 
V do not depend on Pi nor qi — q2, the integration over these variables in Eq. (5.1) leaves 
these V factors intact. This allows us to define a kernel as follows 

lC{P2,t2-ti,z) = ^^iii^^e [ {Pi-P2)S^^\At,Pi,P2,qi-q2), (5.7) 

[Z[L Z)Pq\ Jqi-q2,Pi 

where we have made explicit the independent momentum variables on which S^^^ depends. 
A straightforward calculation of the gaussian integrals in Eq. (5.7) yields 



(5.8) 

Using this definition for /C, one can rewrite Eq. (5.1) as 

= 2g^z(l - z) r dt2 r dti [ V{ka-q2,tL-t2)V{kb-qi+q2,tL-t2) 



X /C(P2, At, z)V{qi -po,ti- to) , (5.9) 



-29- 



where dahard/'i^^po = \ J{Po tPo)\'^ ■ Now, as already explained, the kernel IC{t2 — ti) effec- 
tively restricts the time integrations in Eq. (5.9) to At = t2 — ti < r^^. This is a small time 
interval as compared to the typical values of ti and t2, of order t^, so, in line with our 
previous approximations, we integrate the kernel over At while neglecting the difference 
t2 — ti in the various factors V. We then redefine the kernel after integration over At: 
Eq. (5.8), to be later restored) 



i2—to 



dAt 



1 



p2 
-'2 







2z{l - z)p^ \ cosh^(OAt) 



exp 



2z(l 



ip2 

— ^ tanh(17At) 
- z)Pq^ J 



1 



rt2-to ^ 
z / dAt exp 



iPi 



dAt 



, _^ - tanh(J]At) I 
2z{l-z)p^n ^ '\ 



exp 



2z{i - z)p^n 



tanh(17(t2 - to)) 



1 



(5.10) 



To further simplify the kernel, and again in line with our approximations, we neglect the 
region of integration t2 < r^^^ + to- When t2 — to ^» t^^^, one can use tanh(r2(t2 — to)) ~ 1 
and then the kernel becomes time-independent. After taking the real part in Eq. (5.10) 
and restoring the proper factors from Eq. (5.8), we finally obtain: 



IC{P2,Z,p+ 



z{l - z)p^ 



■Pggiz) sin 





exp 


p2 " 

-^2 






br 




br 



(5.11) 



This kernel generalizes the result obtained in [25] in the eikonal limit. 

Putting everything together, and proceeding to various relabeling (t2 — >• t, -P2 — )• 
p — zq,q2 — 7- p,qi — q), one recovers the expression for the cross-section for quasi- 
instantaneous medium-induced gluon branching given in Eq. (2.21) 

The kernel JC{p — zq, z,Pq) describes the splitting of a gluon with longitudinal momen- 
tum Pq and transverse momentum q into two gluons, one with longitudinal momentum 
fraction z and transverse momentum p, the other with longitudinal momentum fraction 
1 — z and transverse momentum q — p- Note that IC does not depend on p and q separately, 
but only on the combination p — zq, which may be understood as the relative momentum 
of the effective non relativistic two-dimensional motion of the gluons in the transverse 
plane (with q playing the role of the 'center of mass' momentum). We may also write 
p — zq = p'^{vp — Vq) = (p+ — q^){vp-q — Vq), whcrc Vp = p/p'^ etc. are transverse veloci- 
ties. Interestingly, the quantity — is a measure of the actual emission angle in three 
dimensions. To see that, let us introduce the three-dimensional velocities, Vq 
and Vp = {vzp,Vp), with v'^q + Vq = I, etc. Then, {vp — Vq)'^ = 2(1 — cosOp 
similarly for \vq-p — Vq\. Since Eq. (5.11) constrains the value of p — zq to be of order k^^^ 
at most, emission angles are constrained as follows 



(.Vzq,Vq) 

dp, and 



1^ I ~ 
I "pi — 



Vn 



< 



"br 

zq'^ 



— \V, 



q p 



< 



:i 



(5.12) 



These formute make it clear that it is the softest offspring gluon which is emitted at the 
largest angle and hence that controls the geometry of the branching. 
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We may get more precise on the probabilistic interpretation of the kernel /C by inte- 
grating the cross section over suitable elements of phase space. Consider the formula (2.21) 
that we have just derived, 



dVi ~ rtL 



2g'^z{l-z) I dt / V{ka-p,tL-t)V{kb-q + p,tL-t) 



/ * / 



X /C(p - zq, z,p^)V{q -po,t- to) ^l^p^ , (5-13) 
in which it is understood that + — Po ■ From this we may calculate 



, - , d^ka / d^kb (5.14) 

dzdn+ 2 J "J " dkadkbdzdn+ ^ ' 

This quantity represents the cross section for a gluon with initial + momentum in the 
phase space element dfijj" = dp^ / {{2iT)2p^) to split into two gluons, one of which carries 
the fraction zp^ of the initial + momentum. The factor 1/2 is a symmetry factor that 
accounts for the fact that each configuration of 2 identical gluons is counted twice in the 
integration over ka and fcf,- The calculation is easily done and yields 

with 

(5.16) 



dl7+ 7 (27r)2 df]p„ ■ 



From Eq. (5.15), one reads the probability per unit time for a gluon with energy p^ and 
transverse momentum q to produce a splitting with one of the produced gluon carrying a 
momentum fraction zp^ and a transverse momentum p: 

d^=(^^(^-"^'"'^°)- ('-'^^ 
The integration over the transverse momentum yields 

/C(z,P+)= / t\^^^^^^^'') = ^^f^y (5-18) 
J (2vr)^ 27r T^X^.p-^) 

Note that the gluon spectrum produced via a single medium-induced emission may be 
recovered from this formula, by integrating Eq. (5.18) over time up to t/, and multiplying 
by a factor of 2z : 



diV 
'^z 



one emission 



This is recognized as the BDMPSZ spectrum, as expected [8-12]. 

The kernel in Eq. (5.11) may turn negative for large value of the transverse momentum. 
This would not be a serious problem in practice since this occurs for momenta (g > ttAjj^J 
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for which the gaussian in Eq. (5.11) is very small, but it indicates a limitation of the 
approximation that has been used to arrived at Eq. (5.11). In fact, a large transverse 
momentum q S> k^^ cannot be acquired over a time scale ~ t^^ as a result of soft multiple 
scattering. Rather it must be associated with some rare but hard collision, which can be 
more accurately treated in the single scattering approximation — that is, by keeping only 
the term linear in the dipole cross-section in the expansion of the 2-point function (2.14) 
— , but by keeping the logarithmic dependence upon r. When applied to the momen- 
tum broadening probability, this procedure yields for the Fourier transform in Eq. (2.16), 
^(Ap, At) ~ ((7At)/(Ap)^. We expect this procedure to lead to a similar behavior in 1/q^ 
at large q for the kernel IC{q, z,p~^). We leave the detailed discussion of this particular 
point for a subsequent study. 

6 Conclusions 

In this paper, we have provided a complete calculation of medium-induced gluon branching 
in the regime where the gluons that take part in the branching undergo multiple soft 
scattering with the medium. The kernel that describes the branching has been calculated 
as a function of transverse momentum, beyond the eikonal approximation. An important 
conclusion of our study is that the offspring gluons lose color coherence with respect to each 
other on the same time scale as that of the branching process itself. Thus, as soon as they 
are produced, they propagate independently from each other, a picture that holds to within 
corrections of order t^^/L <^ 1. Our explicit proof relies on a large- A'^c approximation but 
we believe that our conclusions remain generally valid. 

In the regime considered in this paper, where the energies of the emitted gluons are 
within the range co^^ <C w <C Wc, a regime dominated by multiple scattering, subsequent 
emissions by these gluons do not interfere with each other. Indeed, the typical duration 
Tj^j. (w) of a branching process is much smaller than the longitudinal extent L of the medium. 
Since the interference effects between the offspring gluons are possible only during this 
branching time t^^^, whereas independent emissions by these gluons can occur anywhere 
along the size L of the medium, the longitudinal phase-space for interference phenomena 
is suppressed compared to the corresponding phase-space for independent emissions by a 
factor T^^{ijj)/L <^ 1. 

The suppression of interference effects is a key ingredient for having independent emis- 
sions. The other key ingredient is that successive emissions do not overlap with each other, 
a situation which, as recalled in the introduction, may occur for sufficiently soft gluons, 
that can be produced at a high rate. Whether a fully probabilistic description of successive 
branchings can be given depends therefore on how well one can control the emission of very 
soft gluons (beyond the apparent infrared divergences produced by the splitting functions, 
and that are easily seen to cancel in the calculation of physical processes). We leave this 
point for a subsequent study. 

There are further limitations of the present calculation that need to be emphasized. 
The calculation that we have presented is valid in a regime where interactions with the 
medium are dominated by soft multiple scattering. While this is a legitimate assumption 
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for the leading particle and the hardest component of its medium-induced radiation, this 
may not be also the case for the very soft gluons at the end of a cascade, especially if their 
energies approach the Bethe-Heitler energy. A proper treatment of this particular region 
would require a more accurate treatment of the single scattering between a gluon and the 
medium, which implies in particular relaxing the eikonal approximation. Also, we have not 
included in our formalism the vacuum-like emissions, that is, the emissions by which the 
energetic gluon that enters the plasma loses its initial (potentially large) virtuality. Such 
emissions are controlled by the standard splitting functions for bremsstrahlung. These are 
generally hard emissions at small angles. In fact, all the emissions at angles much smaller 
than 9c should proceed exactly as in the vacuum with usual destructive interferences leading 
to angular ordering. Hard emissions at larger angles ^ 9c are possible as well, and for 
them interference effects become negligible because of medium effects [21-25]. For hard 
emissions at intermediate angles 9 ^ 9c the situation is more complicated [21, 23], so in 
order to better understand this and also to have a unified description of the in-medium 
jet evolution, it would be very useful to extend our formalism by including vacuum-like 
emissions. 
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A Gluon dynamics in a background field 

In this appendix, we review briefly properties of the gluon propagator in a background 
field A~[x'^,x) which is independent of x~ . To within an inessential term that can be 
ignored, this propagator can be entirely expressed in terms of a 2 + 1 dimensional "scalar" 
propagator (independent of Lorentz indices) , that describes non-relativistic propagation in 
the transverse plane, with playing the role of time. The dependence on Lorentz indices 
is factored out into coefficients that are used to define an effective three-gluon vertex. This 
vertex and the scalar propagator can be used to simplify perturbative calculations for the 
effective non relativistic gluo-dynamics in the transverse plane. 

We work in the light cone gauge = 0, with covariant derivative Z)^ = — igA^, 
and the only non- vanishing component of the background field is A'^ {x^ ,x). In this gauge, 
the free propagator reads, in momentum space, 

G^^(p) = -Go(p)(i^^(p) d^^-'ip)^ 

with n'^ = (n^, n^, n±) = (0, 1, 0), so that n ■ A = A^ = and n ■ p = p^. Note that G^'^ 
is symmetric under the interchange of ^ and z^, and it vanishes if either index is +. The 
non vanishing components are 

Gi^ip) = GoipW', G,\p) = ^Goip), G,-{p) = 2^—G,{p). (A.2) 



n ■ p 



Goip) 



-, (A.l) 
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Figure 11. Representation of the vertex coupling the propagating gluon to the background field A~ . 
The matrix is a matrix of the adjoint representation. The background field A~ is independent 
of x~ , so that there is no transfer of + momentum at the vertex, that is, p'"*" = p"*". 

In the presence of the background field, the propagator is modified. In particular, it is 
no longer diagonal in momentum space (except for the + component of the momentum). 
It is is also a non diagonal matrix in color space. It may be written quite generally as 
G^^{p',p), with the convention that p is incoming and p' is outgoing (and similarly for 
the color indices). That is, one may regard G^^{p',p) as the matrix element (p' b\G^'^\pa), 
which we shall also write {p'\G^^\p) when convenient. To determine the structure of G^'^, 
we write 

ip'b\Gnpa) = ip'\p)SabG^,''ip) + G^'^'(p') (p'6|V,'|pa) <"(p), (A.3) 
with {p'\p) = (2^)4(^(4) (p -p'), and expand T^''' in powers of the background field. We get 

(pmp) = -iipK:''\p)A-ip'-p) 

- iP^K^lp") K{P - P") Go,,'u'ip") {P"K^1P) A-ip" - p) + • • • (A.4) 

Here, {p'\V^^'^\p) stands for V^^^{p' ,p — p', — p), where V is the usual three-gluon vertex, 
defined generally as (with all momenta chosen outgoing), see Fig. II, 

ya7fik,p, q) = -ar"' [g'-'ik - pY + g'^'iP - qT + g'^iq - m ■ (A.5) 

The fact that in Eq. (A.4), or in Eq. (A.3), V''^^'{p',p) is contracted with Go propagators, 
puts constraints on the values of the indices /i and v (see Eqs. (A.2)). A simple analysis 
reveals that the only relevant components of the vertex are 

vi;:{p\p-p',-p) = gf''^g'\p + pY = -'^igp'^in,, g'\ (a.6) 

where we have used p'~^ = p+. It follows that T^'u'{p' ,p) vanishes unless the indices fi and 
V are spatial indices, and furthermore T*-' is diagonal. We shall set 

rij{p',p) = -6ijr{p',p). (A.7) 

Taking this property of T^'^ into account, one can rewrite Eq. (A.3) as 

ip\GZ\p) = ip\p)Sab G'o'ip) - G'o'ip) iPlTM Gtip), (A.8) 
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and, using (A.l), as 

ip\GZ\p) = -ip'\p)Sba G^{p)d''^ - d-'^d^'^ Gob') iPlTbalp)) Goip). (A.9) 
At this point, we define a 'scalar propagator' {p'\Gi,a\p) 

{p\Gba\p) = ip\p)Sba Goip) - Goip') ip'lTbalp) Goip), (A.IO) 
and substitute this into Eq. (A.9) to obtain 

GC(P',P) = ip\p)SbaGoip) [-d^''ip)-d'^\p)d'''ip)]+d'^\p')d'''ip)GUP,p)- (A.ll) 
A direct calculation, using the explicit expression of d^*^ given in Eq. (A.l) reveals that 

ip'lp) Goip) i-d'-'ip) - #Xp)d^"(p)] = -5'^-6''-ip'\p) (A.12) 

We recognize the instantaneous contribution to the gluon propagator in light-cone pertur- 
bation theory. This contact term can be ignored in the present calculation. Thus we are 
left with 

ip'\GZ\p) = d''\p')ip'\Gba\p)d'''ip), (A.13) 

which expresses the gluon propagator as the scalar propagator multiplied by factors that 
carry all information about Lorentz indices. 

The Fourier transform of the scalar propagator, (x|G|y) = G(x,y), is the solution of 
the equation 

- 2igiA- ■ T)d+] Gix, y) = 5(x - y), (A.14) 

an equation that naturally arises when solving the Yang-Mills equations for a fluctuating 
gauge field in the presence of the A~ background [36, 38]. This is easily seen by writing this 
equation as Gq^ + S = 1, with S the self-energy, and then identifying 7" = S — SGqS + • • • . 
Eq. (A.14) can also be written as 

- Vi] Gaix, y) = Six - y)6ab, (A.15) 

where V~ = d~ — igA~ - T, and we used the fact that d~^A~ = 0. Since the background field 
A~ does not depend on x~, the propagator G(x,y) depends on x~ and y~ only through 
the difference x~ — y~ . It is then convenient to introduce a new Green's function 

ix\G'''\y)^ I ^e-^'^^^^--y-^^ix\g''\x+,y+;k+)\y). (A.16) 

A simple calculation reveals that G satisfies the following equation 

V2 



w + 



2k+ 



x\ g'\x+,y+; k+) \y) = iSab^ix^ - y^)Six - y), (A.17) 
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with iV^ = id^ + gA^ . Thus, Q{x^ ,y^; k^) is the propagator of a Schrodinger equation 
for a non relativistic particle of mass k'^ moving in a time dependent potential A~{x'^ ^x), 
with playing the role of the time. It can be written as a path integral 

Pr.e^-/y+ ^*'''[7(x+,y+;r), (A.18) 
with r{y^) = y, r(x^) = x, and U is a Wilson line evaluated along the path r{t) 

C/(x+,y+;r) = Texp|iff^^ di ^"(t, r(i)) T-^j . (A.19) 

From this relation, one easily deduces the following composition property (using the matrix 
notation) 

g{x+, y+; k+) = g{x+, z+; k+) G{z+ , y+- k+). (A.20) 

Consider now the vertex that describes the splitting of a gluon of momentum q into a 
gluon of momentum p and a gluon of momentum q — p- This is given by the general ex- 
pression (A. 5), in which the momentum q is chosen incoming while the other two momenta 
are outgoing: 



vZ\q^ p,q-p) = gr'' [g^'iq + + g'^Q - ^pT - <?'^(2g - pr\ . (A.21) 

A modified vertex is obtained by combining V^^^ with the factors d^^ contained in the 
propagators attached to the vertex. Thus we define 

Cbc = d''i'l)d"'{p)d^''{q-p)V^^l{q,p). (A.22) 

Note that the indices on T are all transverse. Taking into account that the only non 
vanishing contributions of d^* are either of the form d^^ or d^~ , one obtains easily 

r^^^(p - zq; z) = -2gU, |--l-(p _ zqYd^^ + {p - zq^S^^ - -^{p - zg)^5^'| , 

(A.23) 

where z = p'^ /q'^, and we have made explicit the dependence on the single momentum 
p — zq. This momentum has a simple physical interpretation in the effective non relativistic 
dynamics. It can indeed be written as zq~^{vp — Vq), where Vp = p/p'^ and Vq = qjq^ ■ 
That is, p — zq 'xs proportional to the (transverse) velocity of one of the produced gluon 
relative to that of the parent gluon. 

It is convenient also to consider r^^^'^(p — zq\ z) as a matrix connecting single gluon 
states to two gluon states. We set 

(p, 6, i; fc, c, /| r(z) |q, a, i) = {k\q - p)T'^Up " z), (A.24) 

where it is understood that p+ = zg"*", /c+ = {1 — z)q^. Often, we also write ^^aL^P ~ 
zq; z) = gf'^^'^ T^^\p — zq, z), isolating the coupling constant and the color factor from the 
momentum dependent piece. 
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In summary, the dynamics of gluons in the hght cone gauge A'^ = 0, and in the 
presence of a background field that does not depend on x~ , are essentially those of a two 
dimensional non-relativistic field theory. Usual perturbative calculations can be simplified 
by using the propagator Q and the vertex 5/"'"^ Y''^\p — zq, z) that we have discussed in this 
section. Note that the factor l/2p+ which accompanies the propagator Q (see Eq. (A. 16)) 
is not to be included on the external lines. 

B Explicit calculations of n-point functions 

In this section we calculate explicitly the n-point functions that are used in the main text. 
These n-point functions are products of propagators in the amplitude and its complex 
conjugate, i.e. products of ^'s and Q'^\ taken between the same initial time and the same 
final time y"*". At these times, the transverse coordinates in the various propagators may 
take different values. We shall denote generically by X the set of transverse coordinates 
in the various propagators at time x"^, and by Y the set of the corresponding coordinates 
at time . Often, we shall include x"*" and in the sets X and Y , respectively, that is, 
we shall set X = (x"*", X), Y = Y). In case of a single propagator, the same notation 
will be used, with X = {x'^,x) and Y = {y^,y). Further details on the notation will be 
given as we proceed. 

Let us recall that the free propagator reads (we set Go°'{Y,X) = 6^"-Go{Y, X)) 

— )e-^, A(..--.-, (B.1) 

where we have called a; the + component of the momentum, a notation that will be used 
throughout this appendix. We also write ^0(^5 -'^) as the matrix element (^y\QQ{v[^ ^x~)\x) = 
(7/|^o(^^)|^c)- The Fourier transform is 

(q|go(At)|p) = y dyAx^-^^-y ^'^-^ [y\Q^{At)\x) 

= i27rf6{p-q)e-'^^'. (B.2) 

B.l The 2-point function 

The scalar 2-point function S^'^\Y, X\uj) is defined by 

§bb^giab(^j^^Y\uj)g''^{Y,X\uj)) = d""'^ S^'^\Y, X\lu) . (B.3) 

Note that, in order to avoid the proliferation of arguments in the functions that we discuss, 
we are using the same symbols to denote different variables, depending on the 'context', 
that is, depending on the object where these variables appear. Thus, in the propagators 
g{Y,X), X = {x+,x), Y = {y+,y), in g^X,?), X = {x+,x), Y = {y+,y), while in the 
2-point function S^^\Y,X\uj), X = {x^,x,x) and Y = {y~^,y,y). We shall use later a 
similar 'contextual notation' for the 3 and 4-point functions, where X and Y will represent 
collectively larger sets of coordinates (figures, such as Fig. 12 will help to fix any ambiguity 
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Figure 12. Graphical illustration for the 2-point function S''^^^(Y, where X — {x^ ^x,x) and 
Y = {y^,y,y). In parenthesis are given the momenta conjugate to the respective coordinates in 
the (mixed) Fourier representation; the + component of the momenta is denoted here by lu and 
is the same in the amplitude and its complex conjugate. In the present study, a 2-point function 
appears typically as the product of a propagator Q in an amplitude and a propagator in the 
complex conjugate amplitude. These propagators are represented here by oriented lines, with the 
arrow pointing to the right for Q and to the left for 



that may arise). To further simplify the notation, we shall most often omit to indicate 
explicitly the dependence on u. Finally, whenever needed, we shall use a matrix notation, 
keeping separate the dependence on time arguments. Thus we shall occasionally write 
S'(2)(y,X) as the matrix element {Y\S^'^\y+ , x+)\X) with X = {x,x) and Y = {y,y), 
with boldface letters referring to transverse coordinates only. 

By using the path integral representation of the propagators Q and Q"^ one gets, after 
doing the medium average. 



S^^\y,X\l,) = J PrPr exp dt (r^ _ _ ^ J'^ dta{r 



(B.4) 



with paths r{t), r{t) going from r = x,r = x at time to r = y,r = y at time y^, 
respectively. The change of variables 

u = r — f, v = -{r + f), (B.5) 

allows us to rewrite the path integral as 

S^'^\y,X\uj) = J VuVv expjio; dt ii ■ v - dta{u)^. (B.6) 

To evaluate this path integral, we first factorize the dependence on the end points, entirely 
contained in the value of the action for the classical paths. It is easily verified that Euler- 
Lagrange equation for u{t), u = 0, is not modified by the interaction (even though the 



-38- 



individual classical trajectories are no longer straight lines, as in the free case, the difference 
u{t) = r{t) — f[t) remains a linear function of t). This is enough to calculate the classical 
action. We have, in particular, 

/■^^ 1 _ _ _ _ 

J ^ Atu-v= — [y - y - {x - x)][y + y - {x + x)] , (B.7) 

where we have used that it is independent of time and set 

At = y+ - x+. (B.8) 
A standard calculation of the remaining quadratic path integral then yields 

^g„(Y.X)gl{X,Y)o^p{-^ j'^^ dt<T(u(t))l, (B.9) 

where a{u{t)) is evaluated along the classical path. 

By performing the Fourier transform with respect to the transverse coordinates only, 
we obtain the "mixed representation" of the 2-point function, S'^'^\Q , P\uj) where Q = 
{y^,q,q), P = {x~^,p,p). That is, the 'momentum' variables denoted by capital letters 
contains the conjugate to the transverse coordinates, but also the light-cone time. One 
may also write S^'^\Q,P\uj) as the matrix element {Q\S^'^\y~^ , x'^)\P), with P = {p,p), 
Q = (9) q)- We have 

S'^^\Q,P\oj) = y"da;dxdi/dj/e-^'?-^e^P-^e-*P-*e^^-^5(2)(y,X|a;), (B.IO) 

where the signs in the various phase factors follow from our convention for the propagators 
G and that make up S^"^^ . Translational invariance of5(2)(y,X|w) (obvious onEq. (B.9)) 
allows us to perform the integral over the sum of coordinates, leading to a factor (27r)^(5(p — 
p — q + q)- To proceed further, it is convenient to change to relative coordinates. In analogy 
with Eq. (B.5) we define Ux = x — x^ Vx = {x + x)/2, and similarly for y, y. The variable 
conjugate to q—q is Vy—Vx- We are interested in the particular Fourier component for which 
the momenta at are the same in the amplitude and in the complex conjugate amplitude, 
that is q = q. In this case, the variable Vy — Vx drops from the phase factors in Eq. (B.IO), 
and appears only in S^'^\Y, X\uj), where [{y — x)"^ — {y — x)'^] /2 = (vy — Vx){uy — Ux). 
The integration over {vy — Vx) then produces a factor (2Tr At /uj)'^ 6 {uy — Ux), which in turn 
implies that u{t) is independent of time, and also that the 2-point function depends only 
on the difference p — q {p + q being conjugate to Ux — Uy), and not on p and q separately. 
One finally gets 

/AT- 
dUx exp iux -ip-q)-^ (t{Ux) At , 

= (27r) V2)(p - p) V{p - q, At), (B.ll) 
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Figure 13. Graphical illustration for the 3-point function S'^^^ (Y, ^I'i'), with X = {x~^ , Xi,X2, Xq), 
Y = (y"*", yi, y2, yo) f^nd w = {uji,uj2,ujo). The typical 3-point function we deal with in this paper 
is the product of two G (right arrows) and one (left arrow). Also, the gluons 1 and 2 are issued 
from gluon 0, so that conservation of the + momentum implies ljq — ui + uj2- 



where the (5-function expresses momentum conservation after taking into account the condi- 
tion q = q. It is interesting to observe that the requirement q = q has cancelled completely 
the "free" contribution to S^'^\Q,P\u}) (except for the momentum conserving (5-function), 
leaving only in (B.ll) the Fourier transform of the interaction part of S^'^\Y, X\u) (the 
last term in the second line of Eq. (B.9)). This cancellation results in the property of 
5(2)(Q,P|a;) of being a function only of the difference of momenta p — q. In the absence of 
interaction, V{p — q. At) goes over to (27r)^5(p — q) and the dependence on At drops out. 
The quantity 



ViAp, At) = y exp 



lu • Ap At 



(B.12) 



can be interpreted as the probability that a gluon acquires a transverse momentum Ap 
while traversing the medium during a time At. 

B.2 The 3-point function 

The 3-point function S^^\Y,X\uj) that we need is of the form 

fblb2bos{3)^Y,X\u:) = ri'^2ao^gbiai(y^^|^^)gfe2a2^y^^|^^)gtaofeo(X,y|L<Jo)) (B.13) 

where we use the contextual notation introduced earlier, with the variables X and Y having 
different meanings in 5(3) and in Q or in (see the caption of Fig. (13)). By expressing 
the individual propagators in terms of path integrals, and performing the medium average 
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of the Wilson lines, one gets 



^(3) {Y,X\u;) = J Vr2VnVro exp | ^ dt {uirf + cja^i - ojoTq) | 

J ^ dt [a{ri - To) + a{r2 - vq) + a{r2 - ri)] | . 



X exp 



(B.14) 



The integral runs over paths with endpoints rj(x+) = Xi and ri{y^) = Ui {i = 0, 1, 2), and 
the three dipole cross sections correspond to the three possible dipoles that can be formed 
with the three gluons. 

We perform the following change of variables [v can be interpreted as the distance 
between gluon and the "center of mass" of gluons 1 and 2) 

It = n - r2 , V = zri + {I - z)r2 - ro , wi = zujq , a;2 = (1 - z)u!o, 

(B.15) 

so that ri — ro = v — zu and r2 — rQ = v + {l — z)u. For the endpoints, we set = X1—X2, 
Vx = zxi + (1 — z)x2 — xq, and similarly for the y variables. In the variables (B.15), the 
3-point function reads 

r (a,. rv+ 1 



S^^\y, X\u) = j VuVvVro exp | ^ dt [v^ + 2v ■ tq + z{l 

J dt [a{u) + a{v - zu) + a{v + (1 - z)u)] | 



X exp 



(B.16) 



The calculation of the path integrals over ro and v proceeds as for the 2-point function. 
One identifies easily that the classical path v{t) is a straight line, that is, the center of 
mass of gluons 1 and 2 move with respect to gluon at constant velocity. This allows us 
to calculate the kinetic part of the action corresponding to the variables ro and and to 
perform the remaining quadratic path integrals over ro and v. One gets 



2 

exp 



^-.(A. + 2Aro) 



X Jvuexpi^^J dtii^-^^J dt [a{u)+a{v-zu)+a{v + {l-z)u)]^ . 

(B.17) 

where, in the integral, v is to be taken as the classical path v(t), and we have set 

Av = Vy-Vx, Aro = yo-xo, il^o = z{l - z)ujo, (B.18) 

with uiQ the reduced mass for the relative motion of gluons 1 and 2. Note that in the 
absence of interaction with the medium, S^^\Y, X\uj) reduces to the product of three free 
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propagators, as it should. We have indeed 

2 



° ^ ' ' ^ \2TTAtJ 2mAt 



Au 



2 



(27r)3iAt3 



5^(a;iArf+(^2Ar^-a;oAr,-^)^ (B.19) 



with An = Ari — Ar2, and Arj = yi — Xi. The effect of the interaction is, as in the 2 
gluon case discussed above, to produce a damping whenever the size of the dipoles exceed 
l/{qAt). To analyze these effects further, it is convenient to perform a Fourier transform 
with respect to the transverse coordinates. 

The Fourier transform (in the mixed representation) reads 

S'-^^Q , P\uj) = / Q-igi-yi Q-iq2-y2 Q-ipo-xo ^ipi-xi ^ip2-x2 ^iqo-yo S^'^\Y,X\u), 

J XI X2 xo yi y2 yo 

(B.20) 

where translation invariance implies pi + P2 — Qi — (l2 = Po — Qo- It is convenient to 
introduce the following change of variables (with unit Jacobian) 

Xi = {1 - Z)U^ + + Xo, X2 = -ZU^+Vx+Xq, Xq = Xq, (B.21) 

and similarly for the variables yi,y2,yo- The phase factor in Eq. (B.20) becomes then e**^ 
with 

(p = xq{pi +P2 -po) +u^[{l - z)pi - ZP2] +fx(Pl +P2) (B.22) 
- yo{qi + 92 - go) - Uy[{l - z)qi - zq2] - Vy{qi + 92)- (B.23) 

As obvious from its explicit expression (B.16), the 3-point function S^^\Y,X\u) depends 
on Xq and only through the difference xq — yo, which is conjugate to the variable 
[pi +P2 — Po + gi + 92 — go]/2 = 91 + 92 — 90- For the particular Fourier component which 
we are interested in, that for which gi + 92 = go, this variable xq — yo drops from the 
phase factors in Eq. (B.20). One can then perform the corresponding integration, which 
yields a factor {2TxAt/ujQ)'^5{Av). The (5-function 5{Av) = 6{vx — Vy) implies that v{t) is 
independent of time. It also implies that, besides the momentum conserving 5-function, 
S^^\Q,P\u) depends only on the difference po — go rather than on po and go separately. 
We are then left with 

S'^^\Q,P\u) = {2-Kf6^^\pi+p2-po) J du^Auydv^ gm,.(pi-^po)-*^r{'?i-^'7o)+*«.-(Po-qo) 

J Puexpj^ j dt ii^-^^ j dt [a{u)+a{v- zu)+a{v + {l- z)u)\^ , 

(B.24) 

where the (5-function is that of momentum conservation, after taking into account that 

gi + g2 = go- 



-42 - 



To proceed further, we evaluate Eq. (B.24) in the harmonic approximation, Eq. (2.17). 
Then the path integral in the second line of Eq. (B.24) reads simply 



-ZvU^mn+z^Hi-^?] I exp I iuo dt + iojiy^ I , (B.25) 



where we have set 




2 



1 



br 4^2 ' 2 

br 

The calculation of the quadratic integral 



^cff = i:q [(1 - zf + 2^ + 1] 



J{ux,Uy) = y Pit exp I iwo y dt + iw^^it^^ | . (B.26) 
is standard, and yields 

J{u^,Uy) = , ^ ^ ' exp { ] . ' Uul + <) coshQAt - 2u^ ■ uJ } 



47r sinh QAt 1 4 sinh Ai 
47r sinh f^At | 2 sinh OAt 



' . , 2 ^^t^ n2 1 2 n2 

smh (its + itj^j + cosh —^{Us - Uy) 



(B.27) 



with ^ = {1 + i)uj^^ and /c^^ = ^eff^-br- 

The Fourier transform can then be completed. We get 

^fSW^ ^ x2e-f2W , 87r[l + + (1 _ 2)2] _ 2[1+.^ + (1-.)"](po-to)^ 

Q, P = 27r 2^(2) -po ^ ^T-TT — e 3^^* 



A:2^ sinh(r?At) ^] ^ ' MA;2^ coth(J]At/2) ^ ^ 4A;2^ tanh(JlAi/2) j ' 



(B.28) 



where we have set Pi = pi — zpQ and Qi = qi — zqo. Setting S^^\Q,P\uj) = {2tt)'^6^'^\pi+ 
P2—po) S^^\Q , P\uj) , we note that S^^\Q,P\lj) is a function of just three variables, po — 
qo, and the variables -Pi, Qi just defined, and which may be interpreted as the relative 
momenta of the 2-dimension motion at and y+, with po and go playing the role of 
respective 'center of mass' momenta. 

B.3 The path integral needed in the evaluation of the 4-point function 

In this appendix we evaluate the path integral entering the evaluation of the non-factorizable 
piece of the 4-point function (see Eq. (4.18)), which in coordinate space takes the form (see 
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Figure 14. Graphical illustration for the path integral I(Y, X|a;), with X = {x^ ,Xi,X2,Xi,X2), 
Y = (2/^,2/1,2/2,^1,^2) and uj = (a;i,a;2). The 4-point function is the product of two propagators 
Q (for gluons 1 and 2, with + momenta oJi and UJ2, and two propagators corresponding to the 
same gluons in the complex conjugate amplitude. The notation is as in the previous two figures. 

Fig. 14 for notation) 

I{Y,X\u) = y PriPr2PriPr2 exp I ^ y dt (wir-i + W2^^2 - ^1^1 - ^a^^i) | 

J dt [a{ri - fi) + a{r2 - + cr{ri - r2) + a{fi - f2)] | . 

(B.29) 



X exp 



The integral runs over paths with endpoints rj(x+) = Xi and ri{y~^) = yi [i = 1,2), and 
similarly for r,. The factor Ci? in the second exponential of the formula above originates 
from the fact that the average of the Wilson lines that leads to this expression involves 
quark dipoles C\ (see Eq. (4.18)). 

We perform the following change of variables (with unit Jacobian) 



u = ri — r2, u = ri — r2, v 

w 



z{ri - n) + (1 - z){r2 - r2), 
z h (1 - ^)- 



(B.30) 



2 ' ' 2 

with as before uji = zujo,uj2 = (1 — z)uJo, and ljq = z{l — z)a;o- The inverse transformation 
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reads 

,^ , V V 
ri = w + [1 - z)u + -, r2 = w - zu + 

, s _ V _ _ V 

ri = w + [1 - z)u - -, r2 = w-zu--, (B.31) 

We also define 

xi = Wr, + {1 - z)u,j, + -y, X2 = w.j;- zu^ + y, 

*i = + (1 - 2;)ti^ - y, X2 = w^-zux-'^, (B.32) 

and similarly for the y-variables. In the variables (B.30), the path integral reads 
Z{Y,X\u) = j VuVuVvVw exp | "y [{ii^ - ii^) + 2vw]^ 

dt [<T{(l-z){u-u)+v)+a{-z{u-u)+v)+a{u)+a(u)]^ 



X e ' 



(B.33) 



The Euler Lagrange equation for w yields i) = 0, that is, the classical path v{t) is a linear 
function of time. This allows us to calculate the integral over v and and obtain 

Z{Y,X\u) = (^)' e^-oi¥^- j P^.P^ie^/:^*(-^-^^) 



X e ' 



(B.34) 



where Aw = Wy — Wx, Av = Vy — and v{t) = + {Av/At){t — Note that if 

one ignores the second line of Eq. (B.34), one recovers I{Y, X\u) as a product of four free 
propagators, expressed in the variables (B.30). 

At this point, it is convenient to go to the mixed representation, and take a Fourier 
transform with respect to the transverse coordinates. We get 

T{Q, P) = Q-^n-yi Q~ii2-V2 ^iqi yi giq2-y2 ^ipi-xi ^ip2-x2 ^-ipi-xi g-ip2-*2 X{Y, X), 

(B.35) 

where denotes the integration over the 2x8 transverse coordinates, and the variables 

P,Q are those of the mixed representation, e.g., P = {x~^ ,pi,p2,pi,P2)- Translational 
invariance (one can easily verify that I(Y, X) is unchanged in a constant shift of all the 
coordinates) implies P1+P2— Pi— P2 = Qi + Q2 — Qi — Q2- After changing to the variables 
(B.32), and taking into account momentum conservation, one finds that Qi + q2 — Qi — Q2 
is conjugate to — Wy. Now, we are interested in the particular Fourier component 
for which Qi + 92 = Qi + 92- For this component, the dependence of the phase factor in 
Eq. (B.35) on Aw = Wy — drops out. The integration over Aw can then be performed. 
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leading to a factor (27rAt/a;o)^'5(Ai; ), while the integration over Wy + Wx leaves a factor 
(27r)^(^(pi + P2 — Pi — p2)- We are then left with 



i{Q,P)= [ 

Jv 



'VxUxUyUxUy 

X Qtuj:-{{l-z)pi-zp2) ^-iUy{{l-z)qi-zq2) g-Wx-((l-^)pi-zp2) ^iuy{(l-z)qi-zq2) 



X / VuVu e' 2 



X e~^^4^ (B.36) 

In the path integral in the expression above, the endpoints of the paths are Ux,Uy and 
Ux,Uy, respectively, and i; = ■Ua; is a constant. 

We shall evaluate the path integrals over u and u in the harmonic approximation, i.e, 
set Ncna{u) ^ qu^/2. We write the path integral as 



/ 



VuVu e4+ C = Co + A = iy (ii' - (B.37) 



In the harmonic approximation we get 

jCi « --r^q {Aiu^ + u'^)-2{A-l)u-u + 2B{u -u)-v + 2v^} , (B.38) 

with 

A = l + z^ + il- zf, B = l-2z. (B.39) 

The shift of variables 

B B 

u^u - ^ V u^u + i^j—^ V , (B.40) 

where i> is a constant vector, leaves the measure of the path integral unchanged, as well as 
Lq-, and reduces i2i to 

A = {A{v? + «2) _ 2 (A _ x)u ■ u] . (B.41) 

At this point, we perform a boost- like transformation, 

which again leaves £o unchanged, that is, C^iu^u) = Co{u',u'). This leaves also the 
measure of the path integral unchanged. By choosing /3 as a solution of the equation 

/,^ + l-^^ = 0, (B.43) 
one can write Ci in the factorized form 

A = -^97' [^(1 + P') - 2{A - 1) P] {u'^ + u% (B.44) 
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Simple algebra yields 



^^^^ ^ Ail + P')-2{A-1)(3 ^ ^^^^^ ^ ±V2z2 + 2(l-z)2 + l , (B.45) 

where the two signs correspond to the two possible solution of Eq. (B.43). One can then 
write the path integral as 

and the convergence imposes the choice of the positive sign for C{z). This is the product 
of two integral of the type of Eq. (B.26): 

J{u',,u'y)J{u',,u'y) (B.47) 

with here 

(B.48) 

and similarly for J with 



(B.49) 

and the values of u'^, u' u'^, v! can be obtained from the inverse transformation 



1 1^+2(^1)'" 



\U J \-P ^ ) 2(^31)^/ 

At this point, it is clear that the 4-point function that correlates gluons a and h is exponen- 
tially damped on a time scale comparable to that of the 3-point function, more specifically 
as 



C Color algebra 

In this appendix we review the color algebra necessary to justify Eq. (4.6) by replacing 
adjoint Wilson lines by fundamental ones. This manipulation can be made in several 
different ways, here we choose to first replace all adjoint Wilson lines by fundamental ones 
by means of the identity in Eq. (4.5) and then remove all explicit color matrices t's and 
/'s by using the Fierz identity tf^t^; = \6ii6jk — ^^^^ij^ki and the Lie algebra relation 
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= if'^^'^t'^. In order to have shorter expressions we wiU use the shorthand notation 
Ui = U{ri), U{fi) = Ui and similarly for the other coordinates and adjoint Wilson lines. 



2db 



Here we can directly apply the Fierz identity to the pairs of matrices and t''. One can 
easily see that the second term of the identity would remove one of the color matrices 
form the traces, in which case the Wilson lines cancel out and one is left with a trace of a 
fundamental color matrix, which is zero. One therefore obtains 



rUb 



uluit^uluit"" 



Tr 



UlU2t^UlU2f 



(C.2) 



Now all t matrices have different color indices. In order to proceed we need to get rid of 
the / symbols, arriving at 



(C.3) 



Now we have three pairs of t matrices with the same color index and therefore we can make 
use of the Fierz identity three more times. Notice that the contribution from the second 
term of the Fierz identity is always zero since it would eliminate one of the matrices inside 
a commutator. By successively eliminating matrices with the indices d, m, and e, we arrive 
at 

Tr (ulUiUlU2 [f^^f] ulU-2fU\Uit'A - Tr (ulUifUlU2 [f^ ulU-2UlUit 



Tt\U{UiU^U2] Tr I t^f/^t/gt^C/^C/i 1 - Tr I C/JC/it/^C/st" j Tr I C/^C/gt'^C/i 



Tr [U{Uit^U-^U2j Tr (^rC/^C/gt/^C/ij + Tr \U{Uit''U^U2f j Tr \U^U2U{Ui 
Tr (u\UiUlu^ TV (u\Ui^ TV (ulu-^ - TV (ulUiUlU2U\UiUlu^ 



-Tr (ulUiUlU-2U\UiUlU2) + Tr (c/|c/ij Tr (UIU2) Tr (ulU-2U\Ui 
The last line above is exactly the expression in (4.6). 
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